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ABSTOACr 

This  is  an  investigation  of  the  stability  of  self-acting,  gas- lubricated  bearings.  Two 
approaches  to  the  solution  are  presented  and  their  results  are  compared.  Also,  the  relation 
is  discussed  between  the  present  work  and  other,  more  simplified,  methods  available  in  the 
literature.  The  particular  case  of  a  360°  journal  bearing  of  infinite  length  is  treated,  and 
the  chmges  necessary  to  use  the  same  theories  with  other  geometries  are  pointed  out.  Avail¬ 
able  experimental  results  are  collected  and  compared  with  theory. 
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1.  INTRODUCTION 


1.  Imtibilitiei  of  GM-Labric«fd  BaarinM  ind  Rel«tad  Pr<Alew»  -  Reviw*  of  Pertinent  Lit- 
•»tur« 

The  (UvelofBMnt  of  fluid- fila  journal  and  thruat  bearinga  coincided  with  the  advent 
of  the  firat  enginea  in  the  latter  half  of  the  nineteenth  century>  Since  then,  many  contrihu- 
tiona  have  been  made  to  practical  bearing  developeient,  auch  aa  tlioae  of  Towera,  Michell,  and 
Kingsbury  and  to  theoretical  bearing  analysis  by  Petrov,  Reynolds,  Soenerfeld,  end  Harrison. 

In  time,  the  demands  on  the  performance  of  bearings  have  gradually  climbed  to  hi(^  levels, 
especially  during  the  rapid  developamnts  following  lor  Id  War  II.  Recent  progress  in  smdem 
technology  has  stimulated  intensive  research  mtd  development  effort  in  bearings  lubricated, 
not  only  by  conventional  liquid  lubricants  such  as  oil  or  water,  but  also  by  so-called  “exo¬ 
tic”  lubricants  and  gases. 

Apparently  Him  was  first,  in  pointing  out  the  feasibility  of  gaseous  lubrication 
back  in  1854.  Harrison  produced  the  first  theoretical  solutions  for  simple  cases  of  the  govern¬ 
ing  equations  (infinitely  long  slider  snd  journal  bearings).  Gas  bearings  are  now  being  suc¬ 
cessfully  utilized  in  applications  where  conventional  lubricants  and  rolling  element  bearings 
cannot  operate.  Examples  would  be  high  temperature  and  cryogenic  apparatus,  inertial  guidance 
components,  low  friction  sensing  elements  such  as  glaucoma  detectors  and  stress-strain  test¬ 
ing  machines,  nuclear  reactor  circulators,  and  many  others. 

Ihese  applications  take  advantage  of  some  peculiar  properties  of  gas  bearings  such 
as  the  availability  of  the  lubricant  (often  the  ambient  medium  itself),  low  friction  level, 
nearly  total  absence  of  contamination  and  of  apparent  inertia.  However,  together  with  these 
advantages  there  are  many  drawbacks,  such  as  no  boundary  lubrication,  need  for  extremely  ac¬ 
curate  machining  and  for  very  stable  materials,  low  load-carrying  capacity,  a  friction-to- 
load  ratio  higher  than  for  liquids,  susceptibility  to  dust  or  other  impurities  in  the  gas,  and, 
most  serious  of  all,  frequent  and  seemingly  unpredictable  failures  due  to  dynamic  instability 
and  system  resonance. 

According  to  their  principle  of  operation,  gas- lubricated  bearings  can  be  divided 
into  two  groups; 

1.  Externally  pressurized  bearings,  in  which  the  load-carrying  capacity  is  due  to 
the  pressure  of  gas  fed  by  a  source  and  which  then  escapes  through  the  narrow 
slit  of  the  bearing  clearance. 

2.  Self-acting  bearings,  in  which  the  load-carrying  capacity  is  generated  by  the 
relative  motion  of  the  bearing  surfaces. 

The  phase  shift  between  the  line  of  displacement  and  load  direction  combined  with  a 
very  low  level  of  damping  cause  both  these  two  apparently  dissimilar  types  of  bearings  to  be 
susceptible  to  self-excited  instability.  Nanely,  externally  pressurized  bearings  are  prone  to 
the  so-called  “air  hammer’’  phenomenon.  This  is  also  often  encountered  in  valves  and  pressure 
regulators,  and  has  been  the  subject  of  both  experimental  and  theoretical  work  by  Connlet  [7], 
Fisher,  Cherubim,  and  Fuller  [13],  Robinson  and  Sterry  [29],  Allen,  Stokes,  and  Whitley  [1], 
Rothe  [30],  Licht,  Fuller,  and  Stemlicht  [19],  Richardson  [27],  [28],  Licht,  Fuller  nd  Stem- 
licht  [19],  Licht  [21],  and  Licht  and  Elrod  [20]. 

The  present  investigation  is  concerned  with  the  self-excited  instability  which  arises 
in  self-acting  gas- lubricated  journal  bearings,  and  which  is  often  referred  to  as  “half-fre- 
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qumcy  ahirl".  Hm  naaon  for  this  bmo  io  that  this  instability  is  diaractarisad  by  orbital 
mtion  of  tha  shaft  in  tba  diraction  of  its  rotation  and  at  ^roxinataly  half  tha  rotational 
fra<|ioncy  (actually  alaaya  aoawahat  lass  than  half).  This  typa  of  inst^ility  ia  ganerally 
vary  daatnietiva  ainea,  for  valuoa  of  tha  rotational  apaad  highar  than  Uia  thraahold,  the  or¬ 
bital  anplituds  groaa  axtrsnely  rapidly  until  contact  betaaan  shaft  and  bearing  is  achieved. 

Otoe  to  the  poor  boundary  lubrication  properties  of  gases,  instantaneous  failure  is  likely  to 
occur.  In  general  tarns  the  likelihood  of  stability  increases  with  decreasing  load,  speeds, 
and  clearances. 

Experinantal  evidence  on  the  existence  of  this  phenoomion  and  some  of  its  character¬ 
istics  cm  be  found  in  nunerous  publications.  Many  of  these  are  concerned  with  production  and 
development  and  make  an  effort  to  claseify  and  find  remedies  for  instabilities.  In  this  cate¬ 
gory  are  the  work  of  Brix  [3],  Cole  mid  Kerr  [6],  Drescher  [9],  Sixamith  [32],  and  Whitley 
and  Bette  [34].  Other  experimentalists  have  been  interested  in  establishing  a  more  general 
body  of  knowledge  and  have  often  resorted  to  seed-eatpirical  approaches.  In  this  category,  are 
the  valuable  contributions  of  Elwell  [12],  Elwell,  Hooker,  Stemlicbt  [ll],  Fisher,  Chenibim, 
and  Fuller  [13],  and  others. 

Due  to  severe  matheeietical  difficulties,  it  was  not  until  very  recently  that  some 
conclusive  theoretical  work  was  perfoimed.  In  the  case  of  cylindrical  journal  bearings  of 
perfect  geometry  and  infinite  length,  concurrently  with  the  present  work  three  other  approaches 
were  attempted  by  Rentsepis  and  Stemlicht  [26],  Cheng  [S],  mid  Auaman  [2].  Other  valuable 
contributions  in  this  general  field  were  made  by  Gross,  [15],  [17],  Pan  and  Stemlicht  [24], 
and  others.  All  these  authors  found  it  expedient  to  introduce  into  their  analyses  several 
severe  approximations,  the  exact  consequences  of  which  have  not  yet  been  established.  Namely, 
Rentsepis  and  Stemlicht  neglected  the  effect  of  fluid  film  history,  (time  dependent  terms). 
Cheng  used  approximate  Galerkin  analytical  expressions  for  both  steady-state  pressure  and 
trmisient  distributions,  and  Aisman  made  certain  severe  assumptions  on  the  relation  between 
shaft  velocity  components.  Ihese  assisnptions  can  be  justified  only  intuitively,  and  in  a  lim¬ 
ited  manber  of  cases. 

All  above-mentioned  analyses  investigate  stability  with  respect  to  small  perturba¬ 
tions  from  ecyiilibrium.  No  acceptable  approach  has  yet  been  devised  to  solve  the  complete  non¬ 
linear  case  of  finite  oscillations  of  the  shaft  center.  Dtoe  to  the  fact  that  non-linearity  ex¬ 
ists  in  the  space  terms  of  the  governing  equations,  even  the  attainment  of  steady  state  solu¬ 
tions  necessitated  the  use  of  numerical  methods,  such  as  for  the  infinite- length  complete 
journal  bearing  by  Elrod  and  Burgdorfer  [lO],  finite  length  complete  journal  bearings  by 
Raimondi  [25],  flat  and  crowned  sliders  by  Gross  [16],  and  finite  length  partial  cylindrical 
bearings  by  Stevenson  and  Oistelli  [33] . 


2.  lAirpose  and  (Xitline  of  the  Present  Investigation 

Ihe  present  work  is  concerned  with  establiahing  the  ranges  of  parameters  correspond¬ 
ing  to  stable  operation  in  self-acting  gas  bearings.  Ev«i  though  the  treatment  included  is 
limited  to  journal  bearings  of  infinite  length,  the  validity  and  effectiveness  of  the  approadies 
enqiloyed  can  be  extended  without  any  difficulty  in  principle  to  many  other  bearing  geometries. 
Two  methods  of  attack  are  used,  both  of  «4)ich  take  account  of  the  time-dependent  (history) 
terms  in  the  ecgiationa  of  compressible  lubrication,  but  which  differ  in  the  hndling  of  non- 
linearities. 

a)  anill  Ptorturbation  Method; 

Non-linearities  are  eliminated  by  restricting  the  analysis  to  possible  shaft  motions 


-  2  - 


THE  PRANKUN  INSTITUTB  .  M  Bmmtk  Md  Dmfilapmm 


I>A  2049-20 

cnly  iriiila  in  •  vary  narrow  mfa  about  an  aqaiUbriuH  poaition,  Raoontly,  thia  apimach  waa 
auceaaafiilly  aaiployad  by  Lidit  and  Elrod  [20],  in  tha  atudy  of  "  air*h«iaar'*  in  axtamally 
praaauritad  gaa  baaringa.  Tbia  awthod  of  attack  waa  outlinid  in  a  pravioua  raport  [S] . 

b)  Orbit  Program; 

Ihe  conplata  non-linaar  aquationa  partinant  to  tha  caaa  are  intagratad  by  nuaierical 
mathoda  to  obtain  tha  ahaft-cantar  oibita  corraaponding  to  any  apacifiad  aat  of  geoawtrical, 
running,  and  initial  conditiona.  Thia  aolution,  obtained  by  awana  of  a  high  apeed  digital 
computer  program  providea  an  idealited  experimental  rig  which  operatea  obeying  exactly  the  aa- 
aumed  governing  e4piationa.  Ihe  importance  of  the  orbit  program  ia  not  to  be  undereatimated 
aince  detailed  information  from  actual  exparimanta  ia  at  the  preaant  moment  extremely  difficult 
to  obtain.  Ihe  attainment  of  exact  data  ia  impeded  by  interference  from  extraneoua  factora, 
limitationa  on  meaaureawnta,  and  dimutaional  requirementa  at  the  limit  of  preaent  tecimological 
capabilitiea.  Since  it  ia  expenaive  to  eatabliah  complete  atability  mapa  fay  uaing  the  orbit 
progrmn,  ita  main  value  at  preaent  conaiata  of  aerving  aa  a  check  on  the  range  of  applicability 
of  approximate  methoda,  which  then  can  be  uaed  for  the  determination  of  the  thresholda  over  all 
rangea  of  parametera. 
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II.  DEFINITION  OF  THE  PROBLEM  AND  BASIC  EQUATIONS 


1.  Hwwolili*  Ewtioa 

Ihe  following  aquations  apply  to  tha  bahavior  of  the  gaa  in  the  claarmce  space: 
Equation  of  Continuity, 

Equation  of  Momentiw, 

Equation  of  Energy, 

Equation  of  State. 


Aa  it  is  coannn  in  lubrication  theory,  the  ratio  of  the  clearance  to  any  characteristic  linear 
dimension  of  the  bearing  can  always  be  taken  much  smaller  than  unity.  Also  the  pressure  gradi¬ 
ent  is  locally  tangeit  to  the  film  so  that  no  pressure  variations  exist  across  the  clearmce. 
The  effects  of  fluid  inertia  and  gravity  are  neglected  in  conpariaon  to  forces  due  to  viscous 
stresses.  The  fluid  is  taken  to  be  Newtonian  with  constant  viscosity  and  with  a  molecular  men 
free  path  small  with  respect  to  the  clearnce.  Due  to  the  amall  thickness  of  the  film  nd  the 
presence  of  metal  boundaries  it  is  customary  to  asaune  an  isothermal  flow;  this  assiiiq>tion  is 
npported  by  an  order  of  mapiitude  nalysis  presented  by  Elrod  and  Burgdorfer  [10].  The  rotat¬ 
ing  member  is  assumed  to  be  absolutely  rigid  nd  to  possess  a  large  enough  polar  moment  of 
inertia  to  keep  the  rotational  speed  essentially  constnt. 


The  general  lubrication  eqiation,  known  as  “Reynolds’ 
these  assunpticns  in  Appendix  1. 


equation”,  is  derived  from 


_1_ 


_3_ 

W 


(p  l!  .  1  (p  «>  f- 


B(Ph) 

Bt 


B(R») 

B9 


] 


(1) 


Figure  1  represents  the  particular  geometry  treated  in  this  work  and  suggests  the  following 
choices  of  dimensionless  parameters 


c 


c 


hj  c  +  xj  sin  9  +  yj  cos  6 

H  =  —  =  -  =  1  +  X  sin  9  +  Y  cos  6 

c  c 

Then  it  is  natural  to  adopt  the  following  parameters 

A  =  52? 

Pa  c* 


(2a) 

(2b) 

(2c) 

(2d) 

(2e) 
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Definitions  (2n  through 


f)  enable  us  to  write  Reynolds’  equation 


in  the  fonn 


3 

L® 

3 

/  BP  \ 

“  3(m)  Bdn)”! 

— 

(  PH®  “ 

+  *T 

— 

< 

II 

1 

n 

i 

'■  '  "■  f  ■  1 

3? 

V 

R2 

3r) 

\  Br,  / 

-39  ar  J 

for  a  bearing  of  finite  length,  which  reduces  to 


(  V  ^P  \ 

B(PH) 

B(PH)  1 

(  PH®  —  1 

=  A 

+ 

— 

V  B9  / 

BT  J 

for  our  case. 
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(2f) 

(2g) 


(3) 


(4) 


Another  conmcnly  used  fonn  of  Reynolds’  equation  is  obtained  by  adopting 


«/-  =  PH 


(5) 


as  independent  variable;  the  equivalent  of  equation  (4)  is; 


(6) 


Equation  (6)  is  especially  suitable  for  numerical  treatment  due  to  the  fact  that  ^6)  is  a 
much  smoother  function  thmi  Pid),  especially  for  large  values  of  the  bearing  running  para¬ 
meter  A,  and  lower  truncation  errors  are  thus  ensured.  For  the  case  of  isothenwl  films  \p 
physically  represents  the  local  mass  content  of  the  bearing  clearmice. 


Ihe  definition  of  reference  pressure  P,  in  the  case  of  bearings  of  infinite  length 
might  give  rise  to  ambiguities,  and  in  this  work  it  shall  be  handled  in  the  following  way; 
Consider  an  infinite  bearing  as  the  central  region  of  a  very  long  but  finite  bearing  so  that 
the  gas  film  is  in  communication  with  the  ambient.  The  anbient  pressure  can  then  be  chosen  as 
a  meaningful  reference.  Perform  now  a  cyclic  integration  in  6  of  equation  (3).  Die  to  the 
periodic  nature  of  P  and  H  and  with  H  =  H(d,T)  only,  we  have 


ZR^ 


3(PH) 

=  a  ^  — —  de 
ST 


32 

V 


^  p2H2d0  = 


2R2A  3 

l2  3T 


-  S  PHd0 


at  steady  state 

3 

3T  "  ° 


(7) 


(8) 
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■nd 

Ji 

or 

i  =  Ai7  +  B  (W) 

ohere  A  wid  B  are  conatants. 

Since  the  integral  in  (10)  haa  the  aane  value  at  both  enda  of  the  bearing,  it  auat 

be  that 
A  =  0 
and 

B  =  ^  P^H^dd  ia  dependent  of  17.  (11) 

Ihe  nunerical  value  of  B  can  now  be  obtained  at  one  of  the  bearing  enda  of  the  bearing  where, 
by  definition  P  =  1. 

B  =  ^  =  ^(1  +  X  sin  0  +  Y  cos  0]?  =  27r  [  1  + -y  +  Y*)]  (12) 

Equations  (11)  and  (12)  are  called  the  "nass  content  rule”  and  were  firat  introduced  by  Elrod 
and  Burgdorfer  [  10] .  For  bearings  aiiich  are  of  infinite  length  in  the  mathenatical  senae  this 
rule  can  be  applied  also  in  non*8teady  conditions  since  "d/dr)  =  0.  Its  validity  can  be  ex¬ 
tended  to  very  long  but  finite  bearings  by  adopting  the  following  model;  Conaider  the  bear¬ 
ing  is  finite  but  long  enough  to  have  small  77-derivatives  everywhere.  At  sero  speed  the  local 
mass  content  is  everywhere  the  sane;  after  starting,  it  will  vary  slowly  to  its  steady  state 
distribution.  Ihe  bearing  under  study  should  then  be  of  sufficient  length  to  make  the  time 
constant  of  this  process  large  in  comparison  to  the  period  of  any  oscillation  about  the  steady- 
state  position  (see  eq.  (8)).  Use  of  equation  (4)  will  now  ensure  the  absence  of  any  axial 
"leaks"  and  still  retain  meaningful  relation  to  practical  situations. 


S  P*H*d0  =  0 


2.  Equations  of  Motion 


With  the  assumption  of  inflexible  shaft,  constant  ali^iment,  rigid  bearing  mount, 
and  large  polar  moments  of  inertia,  the  equations  of  motion  of  the  shaft  reduce  to  those 
which  apply  to  a  point  mass  (namely  the  shaft  center)  as  affected  by  the  integrated  pressure 
mud  viscous  effects. 

Ihe  fluid  film  exerts  tw>  types  of  forces  on  the  shaft;  nonaal  pressure  and  tangential 
viscous  stress.  Ihe  pressure  force  has  the  following  load  carrying  components  per  unit  axial 
length. 


Pressure 

Force 

F 


x-component  ”  ~  R  ^  Pi  ain  9  S 
y-component  ■  Fy  =  B  ^  Pl  cos  6  d9 


(13) 
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i) 


Hw  frietioul  lore**  4m  to  tlM  viaeotta  ahaar  atraaa  on  tKa  anrfaea  of  tho  rotatinf  ahaft 
can  bo  raaolvad  into  too  porta: 


a)  a  torqao  about  tba  cantor  of  tho  ahaft 

Friction  Torqua  -  t?  f  Sf  ($)  id  (14) 

Ihia  tom  will  not  appaar  in  thia  nalyaio  bocauaa  tho  aaaunption  of  conatant  angular  velocity 
takea  the  place  of  the  third  equation  of  notion  (balance  of  nomenta  about  the  ahaft  center). 

b)  a  reaultant  fmce  through  the  ahaft  center  whidt,  in  turn,  can  be  deccnpoaed 
into  the  coaponenta 


Friction  force 
K 


x-coaq>  *  K,  ~  R  ^  S{  coa  0  id 
y-conp  ■  Ky  =  -  R  (f  Sj  sin  0  id 


(IS) 


Sf  =  M 


3/x 


The  expreaaion  for  the  ahear  atress  ia  (aee  expreaaicn  for  /x  in  Appendis 

^  .  Jt  ±  L  = 

h.  ‘‘I 


2R 


ZR 

30 

on  shaft 

surfacs 

/xflR 

jlI. 

d0  * 

cH 

H  30 

(16) 


Performing  the  integraticna  indicated  in  equetiona  (15)  ae  obtain  (aee  Appendix  2): 
X  ^  P  coa  2^  d&  +  Y  ^  P  sin  20  id)  * 

Y(1  -  Jl  -  e^) 


tiM  int 

!Ll  .  if 

F  2R I  F  F 


PP. 


E*  J  1  - 


(17  a) 


F  2R  I  F  F 


In  ^  A 
3  ^  F 


(X  ^  P  sin  20  d^  +  Y  ^  P  coa  20  d0)  + 
X(1  -  'll  -  €^) 


e®  J  1  -  6^ 


(17b) 


It  ia  evident  that  ell  tenaa  in  aquare  bracketa  but  the  laat  are  of  order  1.  Utiliz¬ 
ing  reaulta  obtained  by  Elrod  and  Burgdorfer  [10]  it  ia  alao  poaaible  to  ahow  that  at  leaat 
for  valuea  of  6  up  to  0.9,  the  laat  tenaa  are  alao  of  order  1.  Therefore  the  net  force  due  to 
friction  ia  of  order  c/R  in  ccepariacn  with  the  normal  preaaure  forcea.  To  retain  conaiatency 
with  the  approximetiona  made  up  to  thia  point,  the  effect  of  the  friction  force  will  be  neglected. 


-  7  . 
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HiMt  tht  •quatioM  of  llBtian  rcdaca  to 


thl 

dt3 

«**yi 


F,  =  R  ^  Pi  •»«  S  i0 

Fy  +  W  =  R  ^  P|CO«  d  d0  *  H 


(18) 

(19) 


where  W  ia  the  external  load  per  unit  axial  length,  and  M  ia  the  rotor  aaaa  per  unit  axial 
length. 


Tranafonaing  equationa  (18)  and  (19)  into  diannaionleaa  f(»m, 


—  =  B  ^  P  ain  «  de 

dT* 


where 


B^Pcoaddd^BL 


(20) 

(21) 


Die  problem  atudied  in  this  report  ia  th«i  defined  to  be  the  determination  of  the 
stability  threshold  for  the  system  described  by  equations  (4)  (or  (6)),  (20)  and  (21). 
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III.  SMALL  PBRTUBtATION  TECINIQUE 


M  they  are  rewritten 

(24,  a) 

(24,  b) 

(24,c) 


and  assigning  numrical  values  to  L  and  A,  we  can,  at  least  in  principle,  find  a  solution  to 
system  (24)  in  the  fora  of  values  for  X  and  Y,  and  a  distribution  P(9);  we  shall  call  this  an 
"equilibrium"  solution  and  denote  its  coordinates  by  a  subscript  zero.  Symbolically, 


A,  L  -*  X„.  Y„,  P,(^). 

(25) 

Physically,  this  process  is  equivalent  to  giving  the  shaft  a  steady  load,  spinning  it  at  a 
fixed  rotational  speed,  and,  by  adding  sufficient  external  damping,  if  necessary,  letting  it 
reach  an  equilibrium  position. 

Let  us  now  consider  small  deviations  from  this  equilibrium  state 
(24)  we  can  use  the  following  expressions; 

!,  SO  that  in  system 

X(T)  =  X„  ♦  x(T) 

(26,  a) 

Y(T)  =  Y„  +  y(T) 

(26,b) 

P(0,T)  =  PJ0)  +  p(d,T) 

(26, c) 

Limiting  ourselves  to  cases  for  which  x,  y,  and  p  and  their  derivatives  are  msall  in  comparison 
to  1,  1,  and  Pq  respectively,  we  undertake  the  problem  of  determining  for  what  range  of  values 
of  the  parameter  B  the  pertwbation  quantities  x,  y,  and  p  have  an  exponentially  decaying  be¬ 
havior  in  time. 

Introchiction  of  definitions  (26)  into  syatem  (24),  use  of  the  above-mentioned  rela¬ 
tion  between  A,  L,  X^,  Y^,  and  P  ,  and  the  neglect  of  powers  and  croasproducts  of  perturbation 
terms  as  quantities  of  higher  orMr,  yields  a  corresponding  set  of  linear  equations 


1.  Linearisation  of  the  Equations 

Is  shall  treat  now  the  syatam  of  equations  (4),  (20),  and  (21) 
below  for  the  purpoae  of  ready  visualisation. 


)  dh 

J  dT* 

=  B  ^  P 

tin 

0  d0 

( 

'  dT* 

=  B  ^  P 

COfl 

0  d0  *  B  L 

Setting 

dh 

BPH 

dT* 

dT* 

BT 

-  0 

-  9  - 
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3»p 


♦  f|  {6)  +  fj  (^)  p  ♦  (5) 


dT 


(27, •) 


dT* 


=  B  ^  p  ain  6  66, 


(27, b) 


d*y  , 

-  ~  B  9  p  coa  5  66, 

dT* 

where 

P«'  Hfl' 

£j  (5)  =  2—  +  3— - 


M,*' 


(27,c) 


(28,  a) 


£3  (5) 


^o'  /PoV 
P  *H  *  '  W  J  ' 


(28, b) 


£4  id) 


PoH„^ 
O  O 


(28, c) 


£5  (0)  = 


A  ain  6 


(28,d) 


A  coa  5 


(28,e) 


£7  (0) 


£  2AP/  3AH,'  I 

\  P,H  ®  H,*  h/  I 


\  ^  I  o 

^ J-  ■  .  -  —  .  —  >  coa  6, 

j  H  *  «•  i 


(28,£) 


-  10  • 


f,  (5) 


THE  PRAMCUN  INSTITUTE  •  LrtaiwiiriM  far  BiiMwfc  and  Dmuhfmnt 

l-A  2049.20 

2AP.*  3H.'P.*  mj 


1  ».  3^.  I 

■  J  —  — —  ♦  — i— _  •  — I  cos  6  * 

\  P,H.»  H.*  H.*  \ 


(2B.g) 


«ther«  a  priae  donotos  difforontiotion  with  respect  to  6. 

Ihe  ayatsai  of  goveminc  equetiens  is  now  linear  atd  the  stability  threshold  can  be 
detensined  by  one  of  the  well  knoan  standard  awthoda  such  as  the  Routh.Hurwits  criterion. 

2.  Datensinstion  of  the  Stability  Threshold 

Taking  the  tiaw  Laplace  transfons  as  defined  by 


6(^,8)  =  /  0(^.T)e-T  dT 


(29) 


of  equations  (27),  we  obtain 


32 


— ^  +  fj  (e)  ^  +  [fj  (9)  +  8  f^  (&)]  p  = 


3^2  -  30 

=  [s  fj  (0)  +  {j  (0)]  ii(s)  +  [s  fg  (0)  +  fg  (0)]  y(s), 

82x  =  B  ^  p  sin  0  d0, 

82y  =  B  ^  p  cos  0  d0. 


(30, a) 

(30, b) 
(30,c) 


The  solution  of  system  (30)  is  subject  to  the  condition  that  p(0,s)  be  periodic  with  period 
27r,  or 


p(0,8)  =  p(0  +  2n,a), 


(31) 


Since  the  coefficients  f^  (0),  i  =  2, . 8  of  equation  (30, a)  are  periodic  with 

period  2n,  condition  (31)  reduces  to 


|i(a,  s)  =  p(a  ♦  a) 


l£ 

30 


0=  a 


=  II 

30 


(32,  a) 
(32,b) 


0  =  a  ♦  ^ 


where  a  is  an  arbitrary  number  which  will  be  chosen  to  be  aero.  Therefore,  the  conditions  to 
be  associated  with  system  (30)  are 


-  11  - 


THE  PRANKUN  INSITTUTE  •  Ubtmorks  for  RmmA  md 


=  p(2>r,s) 


lE 

d0 


is 

id 


l-A  2049-20 
(33, a) 

(33,b) 


(9=0  "  '0  =  2rr 

Siaca  x  nd  y  are  independaat  of  0,  the  aolution  of  equation  (30,  a)  is  of  the  form 

p(9,a)  =  S(0,m)  x(a)  ♦  (X^.a)  y(a).  (34) 


Suhstitution  of  (34)  into  (30,h)  and  (30, c)  will  produce  a  systeai  of  two  equations 
in  X  Mid  y,  the  deteninant  of  the  coefficients  of  which  equated  to  aero  ia  the  so-called 
“characteristic**  equation  of  the  systeai. 


All  X  +  Aij  y  =  0 
Ajj  X  +  Ajj  y  =  0 


(35) 


^1 

Aj2 

A22 

0  “characteristic  eqiation’’ 


(36) 


Ihe  sign  of  the  real  part  of  the  zeroes  of  equation  (36)  indicates  a  decaying  or 
growing  exponential  response  of  the  systen  to  a  perturbation  from  equilibriiaa.  Therefore,  a 
negative  real  part  denotes  stability  and  a  positive  real  part  denotes  instability  and,  in 
principle,  the  problen  is  solved. 


Two  difficulties  arise  in  the  practical  execution  of  this  scheaie;  on  one  hand,  the 

coefficients  f-  (0),  i  =  2 . 8  are  not  known  analytically  but  only  numerically,  and,  on 

the  other  hand,  the  general  form  of  the  solution  of  an  equation  of  the  type  of  (30,  a)  is  not 
known  analytically  but  for  a  few  particular  cases.  Both  of  these  facts  make  it  impossible  to 
obtain  analytical  expressions  for  S{9,»)  and  (X^,s). 

As  a  consequence  the  dependence  of  S  and  Q  on  s  cannot  be  kept  explicit  and  the  prob¬ 
lem  must  be  solved  by  trial  and  error.  In  order  to  simplify  and  reduce  the  volume  of  computa¬ 
tions,  and  since  the  determination  of  the  inatability  threshold,  rather  than  of  the  actual 
response  is  the  object  of  this  analysis,  only  values  of  s  on  the  imaginary  axis  are  used.  We 
may  thus  assume  a  value  of  aero  for  the  real  part  of  **s“  and  determine  by  trial-and-error  the 
value  of  its  imaginary  part  which  will  allow  a  single  value  of  the  dynaadcal  parameter  **B’*  to 
satisfy  both  the  real  and  imaginary  components  of  the  characteristic  equation.  The  following 
procedure  is  adopted: 

a)  Asai^i  a  value  of  s  =  io)  and  solve  equation  (30, a)  with  conditions  (33). 

b)  Obtain  the  characteristic  equation,  and  solve  for  the  values  of  the  parameter  B 
which  are  the  aeroes  of  its  real  and  imaginary  parts.  Qieck  if  the  same  value  of 
B  aatisfiea  both  real  and  imaginary  parts.  If  not,  go  to  a  new  guess  of  s. 
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e)  Wmb  •  MliitioB  is  forndt  sad  B  is  ths  thrsshold  vsl«s  of  the  dynasicsl  psrs* 
■stor,  o)  is  ths  eorrospo^ng  diasnsioaloss  frsqssncy. 

Ihs  sctttsl  seltttion  is  dsrivsd  fnsi  ths  follonac  wHiiiNilstians. 

Substitutisf  exprsssion  (34)  for  p  into  squstioii  (30, s)  and  sopsrating  the  coeffici* 
sots  of  X  and  wa  have 


S*  +  fj  S'  +  (fj  +  s  f^  )  S  =  s  fj  + 

(37. a) 

0*  +  fj  0'  +  (fs  +  a  f^)  0  =  a  fg  + 

(37. b) 

Now  let 

S  =  ifc) 

(38) 

S  =  %  * 

(39, a) 

■nd 

(^  =  Or  +  i  On 

(39,b) 

and  separate  real  from  imaginary  conponents.  System  (37)  then  becomes 

(40.a) 

+  fj  Sr'  +  fj  +  w  f^  %,  =  W  fj 

(40,b) 

1  t  ♦  f,  (^  .  a,  q,  =  fg 

(41.a) 

(41. b) 

Ihe  conditions  on  the  solotion  of  these  equations  are 

,8^(0)  =  Sf^{2^r) 

(42. a) 

Sj,'(0)  =  SR'(2rr) 

(42, b) 

)^(0)  =  Sj,(2rr) 

(42.c) 

iS^'(O)  =  Sj,'(^), 

(42. d) 

j%(0)  =  (^(»r) 

(43.  a) 

|<^,'(0)  =  %’(2n) 

(43. b) 

1(^(0)  =  (^(^) 

f 

(43,c) 

'.(^,'(0)  =  0N’(a>F) 

(43,d) 
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Ssvtrsl  wthedbi  «f  aolution  mr*  tried  for  e^ations  <40  thraii|h  43, d).  A  diacus. 
sion  of  thaae  Mthode  can  ba  carried  out  for  any  ayataa  of  the  type 

i*(x)  +  p(x)a*(x)  *  q(x)x(x)  =  r(x)  (44,  a) 

i(a)  -  i(b)  =  0  (44,b) 

»'(a)  -  t*(b)  =  0  (44,c) 


Iba  aolution  of  ayatau  (44)  ia  unique  bacauaa  the  corraaponding  hoangenaoua  ayataai 
*•  ♦  px*  +  qi  =  0 
x(a)  >  x(b)  =0 
x'(a)  -  x'(b)  =  0 


ia incaapleta"  and  adnita  no  non>trivial  aolution  unleaa  p  or  q  contain  a  free  paruaeter 
•hich  ia  free  to  aaauaw  eiganvaluaa. 


One  uell  knoan  nethod  of  aolution  of  ayaten  (44)  conaiata  of  obtaining  two  independ* 
ent  aolutiona  of  the  hoaxigenaoua  part  of  eq.  (44,a)  by  aolving  the  following  ayaten  numerically; 

(  *1*  ♦  P*i'  ♦  q*i  =  0 


=  0 

=  1 

(45) 

■id 

1  ' 

p*2'  q*2  ==  0 

I  *2^*^ 

=  0 

( 

=  1 

(46) 

Ihen  the 

aolution  of  [44a]  ia; 

X  »2r(x)  * 

xjr(x) 

z  “  Azj 

/  «.)*■’■*  / 

K.)  *■ 

(47) 

where  A,  B  are  conataita 

and  W(x)  ia  the  Vronakiai  of  x^  end  Xj. 

A  and  B  are  detemined  by  uaing  (44,b)  aid  (44,c).  A  atable  "aelf.atarting"  method  auch  as 
that  of  Runge*Kutta.Gill  [31]  cm  be  uaed  for  the  aolution  of  ayatema  (45)  and  (46)  but 
loaaea  of  accuracy  in  the  numerical  proceaa  cm  ariae  due  to  the  fact  that  x,,  and  Xj,  have 
rather  violent  exponential  behavior  aid  aaaune  valuea  which  apai  a  range  of  10^  •  10^.  Nwn 
thia  "exploaion"  happena,  it  ia  very  difficult  to  apply  the  boundary  conditiona  end  retain 
more  that  one  or  two  ai|piificait  fignrea.  Diia  inconvenience  wea  actually  encountered,  aince 
the  problem  waa  initially  prograaad  according  to  thia  acheme.  Attenpta  to  improve  the  accu¬ 
racy  were  fniatrated  eapecielly  in  caaea  involving  large  valuea  of  A  and  e. 
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A  Mcond  Mthod  coiuista  of  aolving  ayatcm  (44)  by  nunerica)  ralaxation  or  by  add* 
ing  a  time  derivative  aueh  aa  eeuld  be  met  in  a  diffuaion  equation  and  nuaerically  develop  an 
ea^totic  aolution.  Thia  nethod  waa  alao  actually  tried  with  very  little  aucceaa.  Dm  be^e 
reaaon  ia  the  folloeing:  Oh  one  hand,  the  problan  haa  ‘'cyclic"  boundary  oonditiona  only, 
therefore,  there  exiata  no  a  priori  knowledge  of  the  general  level  of  tiM  nuBMrieal  valma  of 
the  aolution.  Oh  the  other  hmd,  due  to  the  abaence  of  actual  boundariea  where  the  valuea  of 
dependent  variable  or  ita  derivative  are  fixed,  the  numerical  relaxation  or  diffuaion  proceaa 
poaaeaaea  very  little  internal  daping  ao  that  it  ia  very  difficult  to  control  manerical  ate* 
bility.  Aa  a  conaequmce,  extrenely  anall  "atepa"  are  required  to  prearve  nunerical  ataUlity 
with  conaaquently  vary  long  computation  tinea  to  cIom  the  wide  gap  between  inaccurate  "atart* 
ing  gueaMa”  and  the  final  alution. 

A  third  nethod  waa  then  developed.  It  conaiata  of  perforaing  the  following  atepa: 

a)  Solve  the  two  pr(d>lena 

.  ij*  ♦  p*j'  +  q*j  =  0 

)  *l(a)  =  M 

(  *,(b)  =  M  (48) 

and 

*11*  *  P*Il'  ^  ‘I*!!  ■  >■ 

tll(a)  =  M 

=  M  (49) 

vhere  M  ia  an  arbitrary  number. 

b)  Foim 

*  ~  tjj  *  Czj,  (50) 


and  impose  condition  (44,c) 
*'(a)  =  i'(b) 


to  determine  C. 


We  have 

*Il'(b)  -  *n'(e) 

z  '(a)  -  a  '(b) 


(SI) 


Then,  since  condition  (44,b)  is  satisfied  by  z.,  z..,  and  any  linear  combination  of  the  two 
end  since  z  as  defined  in  equation  (50)  setiafies  (44, a),  we  can  say  that 


z(x) 


zii(x) 


*II 


'(b) 


'(a) 


*11*^*^ 

Zj*(b) 


S|(x) 


(52) 
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i*  the  solution  of  systea  (44).  Qie  to  the  feet  thet  the  solution  of  system  (44)  is  unique 
die  final  result  is  independent  of  the  choice  of  M  as  it  cmi  also  be  easily  ahowi  by  direct 
substitution. 

This  method  has  been  found  to  present  several  advantages  over  the  previous  two.  In 
contrast  with  the  first  method,  it  does  not  let  the  solutions  assume  values  which  are  too  ex* 
traieoua  to  the  problem  because  the  "two  point”  fictitious  boundary  conditions  satisfied  by 
B|  and  i||  more  closely  represent  the  physical  situation  than  those  of  systems  (45)  and  (46). 
Moreover,  it  involves  the  least  amount  of  computations  and  no  numerical  stability  problem. 
Indeed,  problems  of  the  type  of  systems  (48)  or  (49)  cm  be  numerically  solved  in  the  follow* 
ing  way: 


Divide  the  interval  a-b  of  interest  in  N  1  equal  sub-intervals  and  write  the  dif¬ 
ferential  equations  at  each  of  the  N  points  dividing  the  intervals.  Replace  the  dependent 
function  md  its  derivatives  by  finite  difference  approximations;  the  simplest,  md,  moat 
often,  sufficiently  accurate  of  these  approximations  are  obtained  with  the  so-called  “three- 
point  central  difference"  formulae,  according  to  which 


t'(Xj)  ~ 


z''(Xi)  ~ 


z(Xj  +  j)  -  2(Xj  .  j) 

2(b  -  al/N 

z(x-  j)  -  2  z(xj)  +  z(xj  .  j) 
(b  - 


(53) 

(54) 


The  solution  gives  rise  to  N  linear  algebraic  equations  in  N  unknowns 


- 1 

0 

0 

o 

o 

o 

0 

o 

_ 1 

z(xj ) 

■Ri  ■ 

Sjj  822  “23®®*  •  ® 

z(x2) 

«2 

®  “32  “33  “34  ®  •  '  ® 

z(x^) 

= 

•  sseee  a  a 

®  •  •  •  ®  ^-l,N-2  ^-1,N 

^.1 

_  ®  •  •  •  0  0  “n,n-i  “nn 

_Z(X|,()  _ 

Here  the  R.'s  are  known. 


A  matrix  such  as  that  of  the  coefficients  of  equatims  (55)  is  called  "Tridiagonal" 
and  any  algebraic  system  of  N  equations  with  tridiagonal  matrix  cm  be  easily  solved  by  3N 
explicit  equation,  in  one  unknown  each.  It  is  evident,  then,  that  no  numerical  stability  pro¬ 
blem  is  encountered  because  no  iterations  are  necessary,  and  that  no  error  due  to  truncation 
of  an  asymptotic  process  ensues. 

In  the  particular  case  under  consideration,  we  are  not  dealing  with  one  equation, 
but  with  two  systems  of  two  equations  (system  (40)  and  (41),  with  conditions  (42)  md  (43)); 
however,  no  additional  difficulty  in  principle  is  encountered.  Indeed  the  algebraic  equations 
resulting  from  reducing  the  equations  from  differential  to  finite  difference  form  are  of  the 
type 
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(56) 


•i  *i.l  +  *>1  *1  *  “i  *i+l  *  =  *i 

h  Yi.l  *  BiYi  *  K  Vi+l  +  =  ii 

*•»*«  *0'  Vo-  *N+1'  Vn+I* 

Die  recurrence  scheme  which  cm  be  shorn  to  lead  to  the  solution  of  (56)  is  the  fol¬ 


lowing: 


a)  Let  Aj  ~  Bj  Dj  “  Ej  “  0 

and  q  =  x„:  Fj  = 

b)  Use  the  recurrence  relations 


-Yi<=i 


'i+1 


WJi  -  X^z.  ■ 


q+1 


RjYj  •  ZjSj 

WJ.  -  XiZi 


E: 


-w.h. 


i+1 


Zih. 


“i+1 


for 


WiYi 


i 

= 

1  - 

N 

where 

w. 

b. 

+ 

a.  A. 

1 

1 

1  I 

z. 

d. 

+ 

a.R. 

1 

1 

i  i 

R. 

— 

a.C. 

1 

1 

1  i 

X. 

= 

i . 

+ 

f.D- 

1 

1 

1  1 

Yi 

= 

8i 

+ 

fiEi 

Si 

r 

h 

- 

fi^i 

XiCi 


X^Zi  ' 


^i+1 


WiYj  -  X.q 


Fi+1 


*iSi  •  B,Xi 

WiYi  -  XiZi 


c)  Dien  the  desired  soluticns  are  given  by  the  recurrence  relations 
’‘i  \+l  *i+l  ^  yj+j  +  q^.! 


(57) 


(58) 


(59) 


I  Vi  =  °i+l  *i+l  Ej+i  yi+i  +  Fj+i 

i  =  N  -  1 

which  is  started  using  the  knowledge  of  md  Proof  of  the  validity  of  this  solution 
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it  givm  in  Appoidix  3.  Dm  Mthod  of  aolution  fw  aystMia  with  siaple  tridiagonal  ■ntricoa 
CM  ba  obtainad  froa  thia  procedira  aa  a  particular  caaa. 

At  thia  point  tha  aolutiona  for  Sjn  Md  (Md  tharefore  fcur  P)  ara  availabla, 
Md  we  CM  aubatitute  in  aquationa  (30, b)  aM  (30,  c)  to  obtain  the  charactariatic  aquation. 

Define: 


^  cos  d  J|| 


(60,  s) 


{ 


cos  0  Ln  f, 


Die  imaginary  Md  real  part  of  the  characteristic  equation  are: 


(IR^  ■  hh  -  KffJR  ^  (Ir  "  Id)  ^  » 


=  0. 


(60.b) 


(61, a) 


(61, b) 


For  one  particular  case  it  is  possible  to  carry  out  the  solution  malytically.  Dtis 
is  the  case  of  an  unloaded  bearing,  for  which  L  =  0,  Md  correspondingly, 

€„  =  0;  P„(5)  =  1. 

Dien  equations  (40)  Md  (41)  become 
Sa"  -  A  Sh'  +  oj  A  SL  =  A  cos  0 

^  ^  ^  (62) 
-  Al^'  -  wAS^  =  wA  sin  6 


(63) 


-  A(]j^'  +  coAC^  =  -A  sin  9 

(J^*  -  A  (]^'  -  Cl)  A  =  ci>  A  cos  0 

with  boundary  conditions  (42)  Md  (43).  Diese  equations  cm  be  solved  Malytically  Md  yield 

(64) 


(65) 


S|^  =  A  cos  9 

+ 

B  sin  6 

^  =  C  cos  9 

+ 

D  sin  9 

(^  =  B  cos  9 

- 

A  sin  9 

(L  =  D  cos  9 

. 

C  sin  6 
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alMre 

A  =  A  [A?  (c^  -  1)  -  li  /  Dm, 

B  =  (A*  (1  -  J)^  +  (1  ♦  J)]  /  Dm, 

C  =  2w  A*/Dm, 

D  =  wA[A*(l  -  o?)  -  1]/  Dm, 

Dm  =  (1  +  A2)2  +  J2a2  [2(1  -  A*)  +  u?  A^]. 

Evaluation  and  aolution  of  the  dharacteristic  equatim  give  the  following  results: 
a)  root  of  (61,  a) 

D  *  Dm 


(!?• 


B  •  D  +  A  •  C 
b)  roots  of  (61,  b) 


(66) 


■n  -  B  ±  -ior*  -  A*  +  C* 

Bt,  —  =  Den  - 

bJ  ^ 


b2  .  D*  t  a2  .  C* 


(67, a) 
(67, b) 


3.  Presmtation  md  Discussion  of  Results 


Because  of  the  need  for  accurate  steady  state  solutions,  it  was  chosen  to  evaluate 
the  stability  threshold  for  all  the  solutions  presented  by  Elrod  md  Burgdorfer  [10].  The 
pressure  distributions  are  available  as  functions  discretized  at  60  points  and  with  m  accu¬ 
racy  of  six  sipiificmt  figures.  The  Uiivac  I  Computer  CIO  Program  is  still  available  at  The 
Franklin  Institute  if  runs  for  different  parameters  or  more  accuracy  are  required. 

Let  us  define  now 


where  "K"  is  the  value  of  the  product  ffl  at  a  stationary  point  for  the  pressure,  md 
€  =  'Tx^  +  =  eccentricity  ratio 

Any  two  of  the  four  parmaeters  v.  A,  e,  L  ooepletely  define  the  problem.  A  list  of  the  com¬ 
puted  cams  is  shorn  in  Table  1.  For  each  run  a  mt  of  oj's  from  0.0  to  2.0  was  tried  md 
results  were  plotted  as  shorn  by  Figures  8  md  9.  Intermediate  results  such  as  the  functions 
of  S||,  S^,  Q^,  were  printed  out  md  plotted.  An  example  is  shorn  by  Figure  7.  It  cm  be 
sem  that  the  functions  are  cyclic  as  required  by  the  boundary  conditions  thus  confirming  the 
.validity  of  the  adopted  awthods. 
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PifM*  10  eoBtaiM  «  plot  of  th*  thn*  roota  (66),  (67, a),  (67, b).  It  ia  iapoaaible 
to  diatioftiali  tlM  fna  tba  aaaa  obtaiiiad  by  tha  nuawrical  proeedara.  Freai  thia  figure  it  ia 
alae  avUaat  that  tha  aaly  poiata  at  uhich  aaata  aay  of  tha  othar  tao  roota  ia  at  w  =  1 
and  aith  a  thraahold  valua  of  tha  atahility  paraaetar  equal  to  infinity.  Thia  is  than  a  proof 
of  tha  faet  that  at  laaat  in  tha  aaaill,  unloaded  baaringa  of  infinite  length  are  alwaya  un- 
atabla. 


Tha  eanplate  raaulta  ate  praaantad  aa  plotted  in  five  different  manners  in  Figures 
2  through  6  and  listed  in  Table  la.  Figure  2  presents  the  value  of  the  parameter  B  at  the 
threshold  of  instability  abich  is  aaaociatad  with  every  computed  point  on  a  plot  of  = 

(2L)  vs  A  with  linaa  of  constant  v  and  e^.  Figure  3  contains  lines  of  constant  e^,  in  a  plot 
of  A  vs.  =  44/^LB)  and  reprodices  some  of  Cheng’ s  [S]  results.  Figure  5  presents  the 
results  as  plotted  for  a)^*  vs.  C*,  where 


(69) 


(70) 


The  parameters  of  and  C*,  suggested  by  Rentzepis  and  Stemlicht,  are  physically  very  meaning¬ 
ful  because  they  separate  the  effects  of  running  speed  and  clearance.  In  this  figure  a  com¬ 
parison  with  senple  Rentsepis  and  Stemlicht  results  is  made.  Figure  6  presents  a  comparison 
with  some  sample  threshold  values  obtained  by  Auaman.  The  parameters  used  are  A,  £^,  and 


_  M  c 
2  w  R  P, 


2 

TT  B 


(71) 


Strictly  speaking,  conclusions  as  to  on  which  side  of  the  stability  threshold  the  stability 
region  ia  located  cannot  be  reached  by  the  results  of  the  small  perturbation  analysis,  aa  so 
far  described,  but  can  actually  be  obtained  several  ways.  First,  practical  experience  shows 
that  higher  rotational  speeds,  larger  clearmice  to  diameter  ratios,  and  larger  masses,  are  all 
unstabilising  factors,  thus  leading  to  the  conclusion  that  values  of  B  between  zero  mid  the 
threshold  correspond  to  instability.  Second,  it  would  be  possible  to  obtain  a  solution  for  a 
value  of  "s”  close  to  one  of  the  conputed  intersections  of  the  roots  of  the  characteristic 
equation  with  the  imaginary  axis  and  observe  the  sipi  of  the  real  part  of  s  in  crnyunction 
with  the  variation  of  the  value  of  the  stability  parmneter  B.  In  such  a  case  it  should  be  that 
a  decreased  value  of  B  corresponds  to  a  positive  real  part  of  “s".  Third,  the  evaluation  of 
the  finite  response  oibit  by  direct  integration,  which  ia  presented  in  section  V  of  thia  work, 
cmz  provide  the  rniswer  by  evaluating  the  journal  orbit  corresponding  to  values  of  B  on  either 
side  of  the  threshold. 
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IV.  FINITE  SHAFT  ORBIT  EVALUATION 


1.  Finite  Difference  Form  of  Pertinent  Equation. 

Dtis  approach  consiats  of  taking  the  choaen  net  of  governing  equationa,  aaauming 
that  all  functiona  involved  are  “well  behaved”  and  amooth,  and  replacing  the  continuoua 
character  of  the  functiona  over  the  intervala  of  intereat  by  their  valuea  at  diacrete  pointa. 
Ihe  additional  approximation  of  replacing  derivativea  at  each  point  by  algebraic  cand>inationa 
of  the  valuea  of  the  functiona  at  neighboring  pointa  reducea  the  problem  f^om  differential  to 
algebraic.  A  aet  of  integro-differential  equationa  can  be  diacretized  by  meana  of  varioua 
techniques  leading  to  more  or  less  convenient  algebraic  schemes  of  solution.  The  choice  of  a 
technique  is  based  mainly  on  three  considerations;  computation  time,  accuracy,  and  numerical 
stability.  The  last  of  these  should  not  be  confused  with  the  stability  of  the  system  repre¬ 
sented  by  the  equations  in  question,  but  consists  of  a  phenomenon  which  directly  results  from 
the  discretization  of  the  problem.  Therefore,  it  does  not  exist  in  the  differential  system. 
Indeed  the  representation  of  a  function  by  values  at  a  finite  number  of  points  can  only  ■ 
adequately  represent  harmonics  whose  wave-length  is  large  in  comparison  to  the  point  spacing. 
All  higher  harmonics  are  actually  destroyed  and  replaced  by  a  random  distribution  of  round¬ 
off  errors.  As  the  function  is  operated  upon  by  the  system  of  equations  in  question  (for 
example,  as  time  goes  on  in  a  diffusion  equation),  the  influence  of  the  false  higher  harmonics 
leads  to  errors  in  the  lower  components.  If  the  chosen  process  is  not  able  to  dampen  out 
these  errors,  the  true  solution  is  soon  obliterated  and  the  phenomenon  of  numerical  instability 
is  said  to  have  set  in.  Semantically,  the  term  suggests  the  violently  oscillating  and  diverg¬ 
ing  behavior  of  the  values  of  the  dependent  variables  as  observed  in  this  phenomenon.  In  a 
diffusion  equation  numerical  instability  can  be  generally  avoided  by  adopting  sufficiently 
small  time  steps.  Obe  can  also  intuitively  see  how  a  smaller  size  of  discrete  interval  will 
increase  the  accuracy  of  the  approximation.  Therefore,  in  general,  both  the  requirements  of 
stability  and  higher  accuracy  can  be  met  at  the  expense  of  longer  computation  time. 

In  writing  the  approximations  to  the  derivatives  of  the  department  functions,  one 
can  achieve  higher  accuracy  by  involving  more  than  the  point  itself  and  its  two  imnediate 
neighbors.  However,  even  with  5-point  formulae,  great  coiiq>lications  are  usually  encmintered 
in  the  treatment  of  the  boundaries  and  initial  conditions,  as  well  as  because  of  the  presence 
of  an  increased  number  of  unknowns  in  each  equation.  The  treatments  presented  in  this  work 
are  all  concerned  with  three-point  formulae. 


In  the  integration  of  Reynolds'  equation  by  numerical  methods,  either  v/i  =  PH  or 
77  =  was  chosen  as  dependent  variable,  for  two  important  reasons.  First  both  uj  W)  and 
77  (6)  are  smoother  functions  than  P  {&)  and  finite  difference  approximations  of  them  lead  to 
smaller  truncation  errors.  Second  if  P  (0)  is  the  dependent  variable,  Reynolds’  equation  re¬ 
quires  a  knowledge  of  that  is  to  say,  a  knowled^  of  the  motion  of  the  shaft.  This 
requirement  is  inconvenient  since,  in  most  cases,  the  motion  of  the  shaft  is  the  output, 
rather  than  the  input,  to  the  problem.  Therefore,  the  two  forms  of  Reynolds'  equation  of 
interest  are 


3 

r  3-/- 

3H  "I 

“  3i/i  ”1 

— 

- 

-  =A 

— —  +  ——  1 

3(9 

L  39 

39  J 

39  3t  J 

and 

3 

r  H  377 

3H  "1  2  A 

1  377  377  “1 

— 

—  __  . 

77  —  - - 

- + - 

3(9 

L  2  3<9 

39  J  ^/tT 

L39  3t  j 

(72) 


(73) 
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For  the  purpose  ot  general  arguaent  let  ua  call  the  dependent  variable  v  and  divide  the  d  - 
axis  in  intervals  ^  9  and  the  T*axis  in  intervals  A  t.  Ihen  any  of  the  nodes  of  the  space* 
time  grid  can  be  characterised  by  the  pair  of  integers  (i,j)  ehere 


9  =  i  £^e 

and 

T  =  j  A  T 

Take: 

V  (9,  t)  =  7J  l/j-i  ^  +  (1  -  Tj)  l/.j 

(9,  t)  TJ  +  (1  -7))  1  -  Vjj.  j) 

35  2  A5 

- (5,  t)  - - — — — - 

3  5^  te^ 

3  V  V.  j  ^  -  i/.  j 

- (5,  t)  =  - 

3  t  At 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


In  the  typical  case  of  a  progressive  time  integration,  77  and  ^  represent  the  per¬ 
centage  of  participation  of  the  function  at  time  r  =  (j  +  1)  .  A  t  as  compared  with  time 
T  =  j  .  A  T,  It  is  to  be  noticed  that  all  values  of  ^  ^  are  unknown  for  any  i  and  that, 

therefore,  if  77  and  ^  are  different  from  zero  the  finite  difference  equations  replacing 
equations  (72)  and  (73)  contain  three  unknowns  while  only  for 

7]  -  ^  -  o  the  number  of  unknowns  is  reduced  to  one  According  to  the  values  adopted 

for  77  and  ^  the  following  schemes  are  commonly  used: 

a)  77  =  ^  =  o.  Explicit  integration. 

The  finite  difference  problem  is  reduced  in  N  algebraic  equations  in  one  unknown 
each  of  which  can  be  explicitly  solved  for.  Ihis  is  by  far  the  simplest  scheme  from  the 
computational  point  of  view. 

b)  77  =  ^  =  1.  Implicit  integration. 

The  system  now  contains  N  non-linear  algebraic  equations  in  3  unknowns  each.  The 
solution  of  such  a  scheme  can  imply  serious  difficulties. 

c)  T)  ~  ^  Semi -implicit  integration. 

This  method  was  first  used  by  Crank-Nicolson  [8]  and  possesses  very  attractive 
features  which  will  be  discussed  later.  However,  it  also  involves  the  solution  of  a  system 
of  non-linear  algebraic  equations. 

d)  V  =  o,  ^  =  1.  Lee’s  Integration.  (33) 

It  has  some  of  the  advantages  of  implicit  schemes  and  results  in  a  linear  set  of  N 
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algebraic  aquationa  in  three  unknoma  each. 

Other  coabinationa  of  valuea  of  17  and  ^  can  be  need  with  apecial  advantagea  to  be 
derived  in  particular  caaea. 


Aa  far  aa  the  equationa  of  aotion  are  concerned,  the  linearity  sinplifiea  the 
problem  conaiderably. 


d2  X 

=  n  i  P  sin  0  d  0 

Id  T* 

(80) 

1  Y 

■  =  B  4  P  CO®  ^  d  0  +  BL 

(81) 

i<1t2 

Ihe  integrals  on  the  right-hand  aide  can  be  approxiinated  by  the  trapezoidal  rule,  for  it  is 
well  known  that  nothing  is  gained  in  applying  more  complex  quadrature  rules  to  cyclic  inte¬ 
grations  unless  variable  intervals  and  weights  are  used.  Ihe  equations  of  motion  reduce  to: 


A  =  R-A  0  •  A  T  •  i  P.  j  sin  (i  A  0) 

i  =  1  * 

(82,a) 

AX  =  )(  •  AT 

(82, b) 

A  t  =  B  •  A  T  •  A  0  •  i  Pjj  cos  (i  A  0) 

i  =  1  * 

(83, a) 

A Y  =  V  •  AT 

(83,b) 

Ihe  evaluation  of  the  orbits  of  the  shaft  can  then 
running  conditions,  such  as  values  for  A,  L,  B,  taking  a  set 

P.®  (i  =  1 . .  N),  Xj,,  Y  ,  We  use  equations  (82) 

X2,  Y2,  ^2'  ^2'  8°  hack  to  fWynolds*  equation,  and  so  on. 

be  carried  out  by  taking  a  set  of 
of  initial  conditions,  such  as 
and  (83)  and  obtain 

2.  Qioice  of  an  Integration  Method 

A  discussion  of  the  relative  merits  of  methods  a,  b,  c,  d  has  been  published  in 
May  1962  by  Michael  [23]  of  International  Busineas  Machines,  San  Jose,  California.  Pertinent 
considerations  in  the  choice  of  a  method  are  accuracy  and  stability.  Let  us  discuss  ‘:hem 
separately. 

a)  Accuracy  Considerations 

Ihe  truncation  errors  connected  with  the  integration  schemes  presented  above  are 

a)  0  [  (£^e)^  +  AT  ], 

b)  0  [  (A  0)2  +  AT  ], 

c)  0  (  (A  0)2  +  (AT)2  ], 

d)  0  [  (A0)2  +  AT  ]• 

We  see  that  the  error  contributed  by  the  0-discretization  is  the  same  in  all  methods  and  can 
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be  reduced  by  ueitif  •  finer  grid  site.  Ae  for  the  influence  of  the  tine  etep  site  A  T,  we 
eee  that  the  &ank«Nieolson  nethod  ia  better  but  that  thia  advantage  ia  aoon  loat  if  an  in* 
creaaed  AT  ia  uaed,  aa  allowed  by  the  greater  atability  of  thia  method  (aee  auction  on 
"Stability  Conaidarationa”).  Even  thou^  the  order  of  nagnitude  of  the  error  at  each  itera* 
tion  ia  readily  definable,  it  ia  difficult  to  eatimate  the  effect  of  error  propogation. 

Some  phyaical  conaiderationa  might  help  in  connection  with  thia  diacuaaion.  Know* 
ing  that  Reynolda’  equation  (dteya  the  principlea  of  conaervation  of  maaa  and  aince  all  the 
atreamlinea  have  been  aaaumed  to  cloae  on  themaelvea  and  to  be  on  a  plane  normal  to  the  ahaft 
axia.  It  ia  neceaaary  that  the  maaa  content  of  any  given  axial  length  of  bearing  be  con* 
aerved  in  tine.  An  accumilaticn  of  error  would  ahow  an  effect  on  the  maaa  content,  aince  all 
other  diatributiona  are  damped  out  in  the  diffuaion  proceaa.  In  an  infinite  bearing,  there* 
fore,  error  propagation  and  accumulation  can  be  monitored  by  periodic  checking  of  the  integral 

f  PH  A  e  or  f  A  e 

In  addition,  correctiona  can  be  performed  by  multiplying  the  preaaure  diatribution  by  a  factor 
which  makea  the  maaa  content  reaaaune  the  original  value. 

It  ia  very  doubtful  whether  conaiderationa  of  accuracy  alone  will  ever  provide  a 
atrong  enough  argument  in  favor  of  any  one  of  the  propoaed  methoda.  Obvioualy,  thou^,  the 
choice  of  apace  atep  and  time  atep  aiae  will  be  quite  inqmrtant  in  netting  the  general  level 
of  accuracy  of  the  approximation.  It  will  be  shown  later  that  the  adopted  method  provides  an 
accurate  enough  approximation  to  the  exact  solution. 

b)  Numerical  Stability  Considerations 


Ihe  condition  for  numerical  stability  of  a  difference  equation  representing  a  dif¬ 
fusion  process  can  be  obtained  by  several  methods.  In  many  problems,  it  is  very  complicated 
to  find  analytically  a  sufficient  criterion  for  stability  which  gives  threshold  values 
reasonably  close  to  the  actual  ones.  Ihen,  in  order  not  to  waste  computation  tim  by  follow* 
ing  a  safe,  but  very  conservative  criterion,  it  is  convenient  to  proceed  by  trial  and  error. 
Ihe  following  procedure  will  be  followed  to  find  a  threshold  value  of  the  ratio  q  =  ^ 

setting  the  upper  limit  of  the  size  of  A  T  which  will  ensure  numerical  stability  for^^^^ 
a  finite  difference  scheste  using  a  space  interval  equal  to  A  6, 


To  analyze  this  problem,  let  us  consider  the  differential  equation  in  question  as 
a  linearized  system  of  the  type 


3  </f  1 

r  B*'/' 

3*  H  /■ 

3 «/.  \*  3  v/,  3  H 

3V-  “1 

— —  = 

(- 

—  J  -  ""  - 

A - 

BT  A 

L  3 

3  0* 

50/  3030 

30  J 

Define  now 


(85,a) 


(8S,b) 
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(85, c) 


which  are  evaluated  aoaMthere  in  the  interval  in  accordance  with  the  aean  value  theoren.  In 
the  following  F|,  F^,  F3  will  be  regarded  aa  conatant  with  the  juatification  that  the  aolution 
conponenta  cauaing  inatability  vary  at  a  nuch  greater  rate  than  theae  tenaa. 

In  infinite  difference  forai  we  have 


AT 


11 

>1 

fi  [>>  ’ 

^  +  (1  -  17)  Vjj  + 

^2 

r  V 

1 

•-* 

+ 

1 

2  A<9 

L 

F3 

A^ 

-  2  V-. j  ‘ 

(1  -  (V+  1  -  2 

-i'  *  -i'-  i>J- 

(86) 


If  V.J  ia  the  exact  aolution  to  the  aet  of  finite  difference  equationa  (86)  over 
the  apace  and  time  grid,  then  the  actual  aolution  can  be  repreaented  aa  the  aum  of  the  exact 


aolution  and  the  deviationa 

y,  j  ~  V.  j  +  €,  J 

I  I  I 


(87) 


Subatitution  of  (87)  into  (86)  and  uae  of  the  fact  that  aatiafiea  (86)  exactly,  yield  an 
equation  in  the  diaturbance 

* [’ - V-V)  •  * 

+  _i_  r ^  (€ij/i*  -  2  ej  +1  +  «ij_\^)  + 


+  (1  -^) 


Ihe  form  of  (88)  ia  the  aame  aa  that  of  (86)  becauae  of  the  fact  that  (86)  ia 


linear. 


Equation  (88)  can  be  aatiafied  by 


e  (0,  T)  =  X  V  eS'*’  +  '*  "  ‘  K  ^ 

n=l 


(89) 


where  every  term  of  the  aeriea  repreaenta  a  aolution.  Subatituting  equation  (89)  into  (88), 
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1  ^ 


we  have  for  •  repreeenutive  tera: 

(T  +  A  T)  +  'PTb,  5  _  e,  T  +  'R  0  _ 

C  **  0  “ 


Rbfl  , 

»  ^  + 


»  A  ,  i  «,  (T  +  A  T)  +  R  b,  5  T  + 

=  F,  AT  jTj  e^  "  +  (1  -  T))  e  " 

<  r  •„  (T  +  A  T)  +  R  b„  (^  +  A  9)  •„  (T  +  A  T)  +  R  b„  (^  -  A  0)"| 

l’’b  ■*  I 


Fj  at 


2te 

+  (I  -  T>) 

AT 


R  b  +  A  •  T  +  R  b_  ((9  -  A  $) 


+  F, 


A0* 


[•n  T  + 

(  r  (T  +  A  T)  + 

I'L' 


]! 


R“b.  (6*6  6)  .  (T  +  A  T)  +  R  b.  0 


a„  (T  +  A  T)  +  R  b„  (6-  6) 

+  e  "  "  I  + 


,  .  r  •„  T  +  R  b  (^  + 

+  (1  -  ^)  e  "  " 

a.  T  +  R  b  (0  -  A  (9  )"1  I 

Ji 


’1 


A  (9)  „  a.  T  +  R  b  5 

-  2  e  "  "  + 


+  e  " 
or 

a  AT 


1  = 


=  FiAT']77e"  +(l-Tj)V  + 


Fj  at 


2  A  (9 


KBn  AT  + 
' 

[ 


Rb„A0  8„AT-Rb„A0 
n  «  ^  n  n 


+  (1  -  7,)  I  e  ■  " 

AT 

A02 


Rb.  AS  -Rb  Ae 

-  e  "  I  f  + 


']! 


1 


R  b„  A  e  .  _ 

n  _  n  n 

-  2  e  +  e 


^  T  <  r  a  A  T  + 

r  RbAe  „  -RbA^ll 

L'  •“*'  Ji 


AT  a.  iT-’Plb. 


]. 


(90) 


(91) 
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Col lectins  terao  and  uaing  the  coaplex  exponential  definition  of  the  funetiona  ain  (b„  L  6) 
and  coa  (b^  A  9),  «w  have 


-  e*"  •  ~  1  +  7]  Fj  A  T  +  I)  F2  Jq  AT  •  (>1^  ain  b^  A  9)  + 


+  2  ^  F3  q  1  “  coa  b„  A  5 


'll 


=  1  +  (1  -  T))  F.  A  T  +  (1  -  7))  Fj  >17^  ('IT  ain  b,  A  5)  + 


(92) 


e  "  repreaenta  the  ratio  of  error  growth  over  one  atep  A  T  and  ia  dependent  on 
the  pertinent  valuea  of  the  coefficienta  b^.  A  aufficient  condition  for  atability  ia 
obvioualy  that 


+  2  (1  -  f)  F3  q  1  -  coa  b„  A^ 


(93) 


or 

9  b  Ad 

1  -  4  (1  -  #)  F3  q  sin?  -2—  +  0  (A  T) 


•*  I  “  4  ^  F3  q  ain*  -i— —  +  0  (A  T) 


(94) 


or 

»  b  Ad 

1  -  4  (1  -  ^)  F3  q  sin*  -2 - 

2 


1  +  4  #  F3  q  sin* 

Referring  to  equation  (84),  we  see  that 


and  by  the  physical  nature  of  H,  1/),  and  Awe  know  that 
F3  >0 


(95) 


(96) 


(97) 
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Dwn; 

■)  For  ^  =  o  (explicit  integration) 

H  i/»  b  A  9 

-  1  <  1  -  4 -  q  ein*  - -  <  1 

A  2 


or 

H  «//  b.  A  5 

1  -  4 -  q  ein*— —  >  -  1 

A  2 


Ihe  condition  for  numerical  stability  is 


AT  A  A 

(A  5)*  **  2  H  2  P  H* 

b)  and  d)  ^  =  1  (implicit  and  Lee's  Integration) 
1 


-  1  <• 


1  -t-  4  Ft  q  sin 


>  1 


always  satisfied. 

c)  ^  ~  2  Integration) 


„  H  .//  2  X 

1  -  2— “q  sin*  (■*  — ") 
A  2 


-  1  <- 


1  +  2 — — q  sin*  (-0— — ) 


■<  1 


(98) 


(99) 


(100) 


always  satisfied. 


We  conclude  then  that  the  implicit,  Lee’s  and  Crank-Nicolson  integration  SMthods 
are  stable  regardless  of  the  relative  sizes  of  A  T  and  A  6  whereas  the  explicit  Mthod  aaiat 
use  values  of  A  T  limited  by 


AT  < 


2  PH' 


4 


Difficulties  seem  to  arise  from  the  fact  that  this  stability  analysis  gives  a  result  depend¬ 
ent  on  the  value  of  ip  which  is  unknown.  The  obtained  bound  on  A  T,  however,  is  extremely 
useful  in  practice  because  the  upper  bound  of  i/i  is  either  known  by  experience  or  could  be 
calculated  from  the  developing  pressure  profile  and  used  to  adijust  the  adopted  value  of  A  T. 

At  this  point,  it  is  important  to  notice  that  the  apparent  unbounded  stability  of 
methods  b),  c)  and  d)  is  probably  not  actual  because,  for  large  enough  values  of  AT,  most  of 
linearizing  assuiiq>tions  made  in  this  development  cease  to  be  valid.  It  is  definitely  true, 
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hamver,  that  for  tkeaa  Mthoda.  AT  can  be  alloMMl  to  aaauM  valuea  one  or  ttm  ordera  of 
■afnitude  larier  than  thoae  for  the  explicit  aethod.  Even  conaidering  the  ^rp  contraat 
between  the  atability  characteriatic  of  explicit  and  iagtlicit  intention  procedurea,  no 
general  reccanendationa  can  be  aada  on  choice  of  a  aethod  without  conaideration  of  the  par¬ 
ticular  problea  at  hand.  Indeed,  all  inplicit  aethoda  involve  the  aolution  of  ayateaa  of 
simultaneoua  algebraic  equationa,  entailing  long  coaputation  tinea  (partly  overahadowing  the 
advantagea  of  inproved  atability),  nore  ccaq>lex  programa  (with  proportionate  debugging  tinea), 
and  the  poaaibility  of  larger  errora. 

The  aolution  of  aiimtltaneoua  algebraic  equations  can  be  very  simple  in  a  few  par¬ 
ticular  cases,  one  of  which  consists  of  a  set  of  linear  equations  with  a  tridiagonal  matrix, 
mentioned  in  section  IV.  This  matrix  is  actually  encountered  in  the  case  of  a  slider  or  any 
one  dimensional  problem  with  two-point  space  boundary  conditions.  If  the  bearing  is  a  com¬ 
plete  cylinder,  the  matrix  of  the  coefficient  is  tridiagonal  but  for  the  addition  of  the  (I,  N) 
and  (N,  1)  elements  which  already  greatly  complicate  the  solution.  If  the  system  of  algebraic 
equations  is  non-linear,  then  relaxation  methods  must  be  used  with  ensuing  long  computations 
and  possible  truncation  errors. 

For  cases  in  which  the  fluid  film  equation  is  coupled  with  the  dynamical  equations 
of  motion  of  the  bearing  component  parts,  the  problem  of  numerical  stability  must  be  studied 
for  the  whole  system.  Since  any  possible  motion  of  the  parts  is  contained  in  the  variable  \p, 
the  numerical  stability  condition  for  the  finite  difference  form  of  Reynolds’  equation  must 
still  be  satisfied.  It  is  possible,  however,  that  more  stringent  restrictions  must  be  imposed 
on  the  size  of  the  time  step  A  T  so  that  the  number  of  points  that  go  into  describing  an  orbit 
is  sufficiently  large  to  confine  the  error  to  high  order  harmonics.  An  accurate  prediction  of 
this  limiting  value  of  A  T  is  made  difficult  by  the  fact  that,  in  the  equation  of  motion,  all 
points  influence  the  orbit  of  the  center  of  the  shaft  at  the  same  time.  It  will  be  seen  that 
the  response  of  the  shaft  center  to  a  general  set  of  initial  conditions  is  composed  of  two 
main  parts:  a  transient  with  dimensionless  frequencies  often  much  larger  than  unity,  and  a 
smooth  orbit  with  frequency  slightly  below  unity.  In  most  cases  of  small  initial  disturbances, 
the  transient  response  dies  out  and  conclusions  regarding  system  stability  can  be  drawn  from 
the  smooth  orbit.  However,  the  integration  procedure  must  be  stable  while  the  transient  is 
being  felt  if  a  valid  solution  is  desired  for  the  following  orbit.  Since  the  number  of  points 
necessary  to  describe  the  transient  is  generally  rather  large,  it  happens  that  values  of  A  T 
must  be  used  which  are  of  the  same  order  of  magnitude  as  those  imposed  by  numerical  stability 
conditions  for  explicit  integration  methods.  It  is  then  a  waste  to  program  the  lengthy  Lee 
and  Crank-Nicolson  procedures  if  their  chief  advantage  is  not  exploited  over  most  of  the  work. 
The  Crank-Nicolson  method  has  the  added  disadvantage  that  the  error  involved  in  each  step 
increases  as  the  square  of  the  time  step  A  T.  Michael  [23]  indicates  that  this  is  not  a 
problem  in  the  range  of  A  T  characteristic  of  the  situation  at  hand  but  it  is  possible  that 
his  conclusions  were  not  obtained  from  a  general  enough  set  of  tests  to  be  valid  in  all  cases. 

Procedures  b),  c)  and  d)  bring  about  considerable  saving  in  cases  where  the  relative 
motion  of  the  surfaces  is  imposed  as  a  smooth  and  relatively  slowly  varying  function  of  time, 
and  in  the  particular  case  of  no  motion.  Indeed,  in  computational  experiments  designed  to 
assess  the  relative  value  of  the  most  common  integration  procedures,  it  was  found  that  the 
explicit  method  had  approximately  the  same  speed  as  a  relaxation  procedure  for  the  solution  of 
steady  state  pressure  distribution  problems,  whereas  the  Crank-Nicolson  method  was  in  the 
average  50-60  times  faster  in  number  of  steps  and  20-25  times  faster  in  actual  conqtutation 
time.  A  very  adequate  relative  evaluation  of  these  methods  is  presented  by  Michael  [23]. 
Forsythe  and  Wasow  [14]  treat  this  subject  in  general  and  have  some  reservations  about  the 
usefulness  of  the  Crank-Nicolson  method.  Because  of  the  number  of  assumptions  that  must  be 
made  in  evaluating  error  btmds  and  stability  characteristics  of  these  methods  especially  in 
non-linear  cases,  it  is  quite  possible  that  Reynolds'  equation  offers  a  particularly  favorable 
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application  to  aome  method  which  ia  only  of  limited  general  value. 

In  concluaion,  the  explicit  method  wax  used  in  orbit  calculationa  becauae  of  the 
large  number  of  pointa  generally  neceaaary  to  deacribe  the  journal  path.  Ihe  Oank-Nicolson 
method  ia  advised  for  all  calculations  of  steady  pressure  distributions  or  cases  of  smooth 
enough  motion  (by  being  externally  imposed  or  for  cases  of  Urge  external  damping). 


3.  Description  of  Special  Features  of  Numerical  Procedures. 

In  comparison  to  the  computation  time  necessary  to  perform  the  iteration  of 
Reynolds’  equation,  even  when  as  few  as  30  circumferential  points  are  used,  the  dynamic 
equations  are  treated  in  zero  time  for  all  practical  purposes.  Many  complex  features  can 
therefore  be  incorporated  with  an  insignificant  sacrifice  in  computer  time. 

One  obvious  generalization  is  the  consideration  of  an  arbitrary  rotating  unbalance 
force.  It  has  been  noticed  by  some  experimental  investigators  that,  in  the  case  of  gas 
bearings  working  with  light  external  loadings,  a  limited  amount  of  unbalance  prevents  failure 
by  ha  If- frequency  whirl.  Unbalance  forces,  of  course,  can  always  be  made  large  enough  to 
cause  failure.  Provision  for  the  unbalance  feature  is  incorporated  in  the  program. 

A  second  special  feature  was  incorporated  for  the  sole  purpose  of  shortening  compu¬ 
tation  time  in  evaluations  of  the  stability  threshold.  When  the  value  of  the  stability 
parameter  B  is  close  to  the  instability  threshold  the  rate  of  convergence  or  divergence  of 
the  journal  center  orbit  is  small.  If  the  orbit  is  now  perturbed  by  transient  oscillations  of 
higher  frequency,  it  becomes  difficult  to  estimate  whether  the  case  is  stable  or  unstable 
before  a  large  number  of  orbits  have  been  accumulated.  Since  most  conclusions  are  drawn  from 
the  behavior  of  the  shaft  after  the  high  frequency  transient  subsides,  a  discriminating 
artificial  damper  is  applied  to  the  shaft.  This  device  operates  as  a  linear  viscous  damper 
and  acts  with  a  force  which  is  proportional  to  the  vectorial  difference  between  the  actual 
journal  center  velocity  and  the  velocity  corresponding  to  half- frequency  circular  or  elliptic 
whirl.  The  artificial  damper  is  applied  at  the  beginning  of  an  orbit  and  is  taken  off  after  a 
specifiable  number  of  steps  as  dictated  by  experience.  The  purpose  is  to  smooth  out  the  high 
frequency  transient  response  in  as  small  a  number  of  steps  as  possible.  It  could  be  inter¬ 
preted  as  letting  the  computer  find  by  itself  a  set  of  “smooth"  initial  conditions. 

The  orbit  programs  are  all  versatile  enough  to  accept  changes  of  physical  and 
running  parameters  while  in  orbit,  thus  being  able  to  follow  an  acceleration  procedure  such 
as  at  “start-up”  or  a  load  variation  in  time. 

Numerical  instability  detectors  and  automatic  time  step  changes  are  also  incorpor¬ 
ated. 


Some  of  the  outputs  are  included  in  this  paper.  The  pressure  distributions  are  ob¬ 
tained  only  at  intervals  of  time  which  are  specified  by  the  input.  This  selection  is  made  in 
order  to  eliminate  wasteful  output  time,  but  the  possibility  of  following  the  pressure  history 
is  left.  This  feature  is  particularly  important  in  connection  with  the  assessment  of  the 
relative  value  and  range  of  application  of  approximate  theories.  For  example,  the  question 
of  the  relative  size  of  the  term  H  with  respect  to  P  and  other  terms  in  Reynolds’ 
equation  can  be  settled  by  direct  computation  from  the  values  of  the  pressure  at  every  point 
of  the  0,  T  grid. 
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4.  vA  Ornmylipp  9i  RffitvUti 

The  whirl  problem  in  cylincricel  journal  bearings  was  attacked  in  steps  in  order  to 
acquire  familiarity  with  available  numerical  techniques  and  develop  and  debug  any  new  methods 
that  might  become  necessary.  Ihe  first  problem  solved  was  the  case  of  the  flat  slider  of 
infinite  lateral  extent  with  compressible  lubrication.  Ihe  geometry  under  consideration  is 
shown  in  Fig.  11.  Ihis  problem  can  be  solved  exactly  by  analytical  methods  and  can  provide 
an  instructive  check  on  the  accuracy  of  the  numerical  techniques. 

For  this  case,  at  steady  state,  Reynolds’  equation  is 

d  /  d  P\  d  (PH) 

—  (PH^  -  )=  A  - 

d0\  i6  /  d  e  (101) 

with  6  -  ^  ^  o-W  (straight  line)  and  can  be  solved  by  the  following  procedure: 
let  0  =  P  H  and  integrate  once 

d  (4>m 

^  =  A 1^  +  const .  ( 102 ) 

d  e 

or,  since  d  0  =  a  d  H 

Aa(/.  +  i//2+(.  =  (103) 

Cl  n 

Ihis  equation  can  be  separated  as 

\p  d  \p  d  H 

C  +  Aai/^  +  i//^  H  (104) 

Now  let  -  A  a  =  9j  +  q^  , 

and  C  =  •  qj 

Then: 


ip  1  ”  </)  'Z'  ” 

+  AaiZi  +  i//^  )_</<“  'p  ~  ^2  — 


Integration  of  (104)  yields: 
9. 


>Z-9i| 


’1 


=  E  H  '^1 


(9,  -  9^) 


'p-% 


where  E,  9^,  and  92  can  be  found  from  the  conditions 


(105) 


(106) 
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P  leading  edge  =  1 
P  trailing  edge  =  1 


=  -  A  a 


(107) 


Figures  12  and  13  show  results  of  computations  for  two  slider  bearing  cases.  Ihese  calcula* 
tions  were  made  first  by  desk  calculator,  and  later  with  the  help  of  an  IBM  650  computer.  Ihe 
initial  pressure  distribution  was  taken  to  be  ambient  throughout  to  that  this  situation 
resembles  the  problem  of  a  suddenly  accelerated  wall.  Besides  the  steady-state  results  and 
the  exact  solution,  Figures  12  and  13  present  pressure  distribution  curves  at  intermediate 
values  of  time  corresponding  to  explicit  integration.  It  should  be  noticed  that  high  accuracy 
can  be  obtained  by  this  procedure,  but  at  the  expense  of  a  very  large  number  of  steps  (in  the 
neighborhood  of  150-200  to  reduce  the  truncation  error  to  a  fraction  of  10*^).  Ihe  Cirank- 
Nicolson  method  applied  to  the  same  cases  achieves  the  same  results  in  approximately  10  steps. 
Figure  12  also  shows  results  obtained  by  Gross  [16]  by  relaxation  methods.  This  plot  was  drawn 
taking  points  from  graphical  data,  so  that  the  apparent  slight  discrepancy  between  Gross’  re- 
suits  and  the  exact  solution  might  be  due  to  reading  errors,  rather  than  high  truncation 
values.  It  is  remarkable,  though,  that  the  computing  times  quoted  by  Gross  as  necessary  for 
relaxation  solutions  are  commensurate  with  explicit  integration  times. 


After  the  successful  completion  of  many  slider-bearing  runs,  the  complete  journal- 
bearing  dynamic  problem  was  programmed  for  a  Univac  I  computer  at  the  Franklin  Institute.  The 
first  successful  set  of  orbits  were  obtained  with  programs  employing  P,  Rl,  and  P^^as  inde¬ 
pendent  variables.  It  was  learned  that  the  use  of  P  was  definitely  detrimental  to  the  accuracy 
of  the  solution,  especially  when  high  values  of  the  running  parameter  A  were  involved.  With  an 
average  speed  of  5  seconds  per  time  step  it  soon  became  a  problem  to  keep  computer  time  down 
to  reasonable  levels.  Indeed,  it  had  been  hoped  that  fifteen  time  steps  would  suffice  to  de¬ 
scribe  one  ha If- frequency  orbit  since  at  the  time  there  was  no  published  evidence  on  the  exist¬ 
ence  of  high  frequency  transients.  It  became  then  necessary  to  evaluate  approximately  200 
points  per  ha If- frequency  orbit.  Moreover,  because  of  the  very  limited  fast  memory  capacity  of 
the  Univac  I  computer,  no  more  sophisticated  integration  methods  could  be  used. 

Fortunately  then,  the  Bureau  of  Ships,  U.  S.  Navy  Department  made  available  some 
time  on  their  Remington  Rand  LARC  and  191  7090  high-speed  computers  at  the  David  Taylor  Model 
Basin  in  Washington.  The  LARC.  was  chosen  first,  but  the  program  had  to  be  written  in  machine 
language,  since  no  working  pseudo-code  was  available  at  the  time.  Consequent  debugging  diffi¬ 
culties  discouraged  further  use  of  that  machine. 


All  subsequent  work  was  conducted  on  the  194  7090.  The  fORTRAN  pseudo-code  connected 
with  the  very  efficient  Bell  systcsm  BE-Sys-3  and  the  very  helpful  staff  of  the  Applied  Mathe¬ 
matics  Laboratory  of  D.T.M.B.  made  programming,  debugging,  and  running  very  simple  and  expedi¬ 
tious.  The  remarkably  high  internal  speed  of  this  modern  computer  (access  time  of  2.4  micro¬ 
seconds,  floating  add  time  14.4  microseconds)  made  possible  an  increase  in  integration  frequency 
to  1600-2000  steps/minute  for  a  30  point  space  grid.  Thus,  seven  or  eight  half- frequency  orbits 
can  be  obtained  in  one  minute;  this  is  generally  sufficient  to  decide  the  stability  of  a  case. 

A  very  useful  feature  of  the  D.T.M.B.  facility  is  the  availability  of  a  General 
Dynamics  CHARACTRQN  which  can  be  programned  automatically  by  the  IBM  7090  to  plot  any  array  of 
X-Y  pairs  on  a  cathode-ray  tube  and  issue  a  microfilm  photograph  of  it.  This  eliminated  most 
of  the  cumbersome  hsnd  plotting  that  was  previously  necessary. 

The  first  check  on  the  7090  program  was  provided  by  a  duplication  of  Univac  results. 
The  only  existing  differences  consisted  of  machine  round-off  errors  which  are  due  to  the  fact 
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that  the  acbwl  7090  word  length  ia  one  digit  shorten,  and  that  the  Uiivac  progrm  waa  written 
in  fixed  point. 

Ohe  of  the  interesting  results  obtained  with  the  shaft  position  held  fixed  was  the 
slight  overshoot  of  the  pressure  distribution  over  and  under  the  steady  state  pressure  profile. 
Ihis  behavior  is  shown  in  Figure  14  and  seems  to  be  dependent  on  the  nuining  parameter  A. 

Namely  high  values  of  A  will  correspond  to  more  ample  and  prolonged  oscillations.  Figure  14 
represents  pressure  profiles  at  equally  spaced  time  intervals.  Since  each  time  interval  cor¬ 
responds  to  100  steps,  it  can  be  seen  once  again  that  the  explicit  method  is  rather  inefficient 
in  obtaining  steady  state  results. 

One  of  the  most  important  tests  to  perform  on  the  finite-orbit  technique  is  to  ex¬ 
amine  if  two  different  values  of  the  time-step  size  AT,  both  such  as  to  ensure  numerical  sta¬ 
bility,  produce  the  same  orbit  under  the  same  initial  and  running  conditions.  This  is  clearly 
shown  in  tabulations  and  orbit  plots  of  Figures  15  through  20.  Figure  19  contains  100  X-Y 
purs  for  a  bearing  running  under  the  conditions  of  Ihin  13  and  with 

B  =  .45 

AT  =  .026 

started  from  the  steady  state  position  of  X^  ”  .1563,  =  .57928  and  with 

=  -01 

Yo  =  .01 

The  listing  of  every  step  on  Figure  19  corresponds  to  every  second  step  on  Figure 
20  which  is  for  the  same  case  with  AT  =  .013.  It  can  be  seen  that  agreement  to  approximately 
five  figures  is  attained.  Figures  17  and  18  show  the  same  results  for  the  pressure  profiles, 
although  data  for  exactly  double  the  number  of  steps  are  not  available.  Other  pieces  of  evi¬ 
dence  to  the  same  extent  are  available  for  examination  in  the  A2049  file  of  The  Franklin  Insti- 
tu  te. 


Another  importmit  verification  of  the  validity  of  the  finite  orbit  program  is  furn¬ 
ished  by  the  agreement  between  steady-state  positions  and  pressure  profiles  obtained  by  Elrod 
and  Burgdorfer  [10]  and  the  values  resulting  from  the  settling  of  the  shaft  center  at  an  equ¬ 
ilibrium  point  in  stable  cases.  The  validity  of  this  argument  is  particularly  evident  when 
one  notices  that  EUrod  and  Birgdorfer  obtained  their  results  by  imposing  the  geometrical  con¬ 
figuration  and  obtained  the  load  magnitude.  A,  and  mass  content  as  results,  whereas  the  orbit 
program  imposes  the  load.  A,  and  the  mass  content  and  finds  a  geometrical  configuration  cor¬ 
responding  to  them.  Thus  the  two  methods  are  completely  independent. 

Choosing  as  a  sample  comparison  Elrod  and  Burgdorfer’ s  case  of 


V  -  0. 5 
e  =  0.6 

and  the  corresponding  orbit  progrmn  for 
A  =  1.460 
L  =  1.9934 
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we  obtain  the  values  of  Table  II  for  i/*  =  M  distributions  at  equilibrium.  Elrod  and  Burg* 
dorfer  use  sets  of  points  which  are  evenly  ^^aced  in  the  distorted  spatial  variable  /8  related 
to  d  by 


cos  6 


€  -  coa  /3 
e  cos  /3  -  I 


(108) 


Iherefore,  a  direct  nusterical  comparison  of  the  two  distributims  is  not  possible 
without  interpolation.  Figure  21  contains  the  two  plots  of  those  distributions  and  shows 
conclusively  excellent  agreement.  The  orbit  which  approaches  this  equilibrium  point  is  shown 
in  Figure  22.  It  should  be  pointed  out  now  that  in  all  CHARACTRON-plotted  orbits  only  one  out 
of  every  four  X*Y  pairs  is  registered  and  a  straight  line  drawn  between  consecutive  points. 

As  a  result  some  orbits  appear  to  be  less  smooth  than  they  are  actually  prochiced  by  the  com¬ 
puter  integration. 

The  foregoing  proofs  of  accuracy  of  the  numerical  procedure  are  deemed  to  have  suf¬ 
ficiently  established  the  validity  of  the  results  «hich  are  hereinafter  discussed. 

C^ceming  the  general  behavior  of  the  solutions  -  it  was  previously  mentioned  that, 
for  values  of  the  stability  parameter  above  the  critical  (stable),  the  response  is  composed 
of  two  simultaneous  parts:  a  rather  smooth  precession  in  the  same  direction  as  the  rotation 
of  the  journal  at  a  frequency  slightly  below  D/2  and  a  higher  frequency  component.  The  fre¬ 
quency  of  the  overriding  component  is  directly  related  to  the  value  of  B,  Indeed,  there  is 
evidence  to  show  that,  at  least  for  cases  far  from  the  instability  threshold,  the  overriding 
frequency  increases  proportionately  to  the  square  root  of  B.  Instability  sets  in  «hen  the 
value  of  B  is  low  enough  to  give  rise  to  overriding  components,  the  frequency  of  which  is 
commensurate  with  rV2.  It  would  that  appear  as  though  the  actual  bearing  response  is  the  over¬ 
riding  oscillation  with  its  frequency  and  damping  characteristics  directly  controlled  by  the 
film  paraneters,  while  the  half- frequency  component  is  a  pseudo-natural  frequency  to  which 
the  bearing  cannot  react.  Then,  when  the  bearing- response  frequency  is  of  the  order  of  D/2, 
the  resonance  is  excited  and  divergence  occurs.  For  unstable  cases,  the  rate  of  divergence  is 
also  controlled  by  the  value  of  B  and  the  half  frequency  component  disappears.  Therefore,  fre¬ 
quencies  much  below  D/2  are  possible  in  a  bearing,  although  it  is  very  difficult  to  obtain 
experimental  evidence  of  them  -  probably  due  to  the  fact  that  the  rate  of  divergence  is  so 
rapid  that  failure  would  always  ensue. 

If  this  model  is  accepted  for  the  phenomenon  of  self-excited  instability,  a  parallel 
can  be  drawn  with  the  commonly  known  case  of  "oil  whip”.  This  tenn  refers  to  the  restmance  of 
the  elastic  response  frequency  of  a  shaft  with  the  half  frequency  component  of  conplete  oil 
journal  bearings. 

This  interpretation  would  then  attribute  the  overriding  response  to  the  "squeeze- 
film"  effect  due  to  the  trapping  of  gas  between  the  moving  shaft  and  the  bearing.  It  can  be 
easily  visualized  that  trapped  gas  gives  rise  to  an  elastic  restoring  force  which,  coupled 
with  the  mass  of  journal  through  the  parameter  B>  ho*  a  characteristic  response  frequency. 
Trapping,  however,  is  not  complete  and  escaping  gas  produced  the  viscous  shear  dmping  that 
slowly  eliminates  the  overriding  response.  In  long  journal  bearings  gas  can  escape  only  in 
the  circumferential  directim  whereas  more  and  more  axial  flow  is  possible  as  the  bearing 
length  decreases.  This  justifies  the  high  level  of  dmping  of  short  journal  bearings. 

The  effect  of  load  can  also  be  justified  by  means  of  this  interpretation.  Indeed, 
keeping  A  constant  and  increasing  the  load  L  causes  an  increase  of  eccaitricity  ratio,  thus 
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ma^iifying  the  elastic  to  deaping  ratio  of  the  squeeie-filn  effect.  Thua.  for  any  given  value 
of  B.  the  filai  response  frequmcy  is  hi|^er  and  the  threshold  of  instability  occurs  at  a  lower 
value  of  B*  I"  fact,  load  is  seen  to  have  a  stabilising  influence,  all  other  paraneters  being 
kept  the  sane. 

The  daBy>ing  characteristics  of  the  film  show  a  dependence  (Hi  the  bearing  parameter 
A.  Namely,  high  values  of  A  give  rise  to  smaller  danping  with  cansec|uent  prolonged  duration 
of  the  trmiaient  responses. 

Many  other  bearing  response  characteristics  might  be  derived  from  examination  of 
the  orbit  program  results.  Much  material  is  already  available  at  The  Franklin  Institute  for 
anyone  interested  in  further  studies,  and  plans  have  been  made  to  exploit  extensively  the 
orbit  program  within  the  sane  O.N.R.  contract  that  is  responsible  for  its  creation.  Here  we 
shall  discuss  results  of  the  finite  orbit  program  vdiich  are  pertinent  to  the  definition  of 
the  stability  region  of  infinite  gas  bearing,  so  that  a  check  may  be  provided  to  the  conclu* 
sicxis  derived  from  the  small  perturbation  approach  of  Section  IV* 

A  set  of  nine  cases  was  run  corresponding  to  iVms  7,  8,  12,  13,  16,  17,  18,  19,  20. 
These  are  in  excellent  agreement  with  the  small  perturbation  results.  Perhaps  the  term  "excel¬ 
lent  agreement”  should  be  qualified.  Indeed,  several  factors  contribute  to  some  uncertainty  in 
the  determination  of  the  threshold  of  instability  by  means  of  the  finite  orbit  program.  First, 
for  values  of  B  close  to  the  threshold  the  rate  of  convergence  or  divergence  is  so  low  that 
increasingly  higher  nusbers  of  orbits  are  necessary  to  interpret  each  case  as  the  threshold  is 
second,  the  threshold  value  of  B  is  not  unique  since  the  size  of  the  initial  disturbance  has 
an  effect  on  the  result.  Indeed,  the  case  studied  to-date  show  once  again  that  perturbation 
analyses  of  stability  problems  are  only  able  to  produce  necessary,  but  not  sufficient,  criteria 
for  stability.  For  certain  rrniges  of  values  of  B  the  system  is  stable  with  respect  to  small 
disturbances  while  it  is  unstable  for  large  ones.  This  phenomenon  is  widely  known  «id  has  often 
been  proven  experimentally.  Licht’s  test  apparatus,  still  operational  at  The  Franklin  Institute 
[22],  very  clearly  demonstrates  that  non-linear  effects  are  often  sufficient  to  overpower  linear 
ones  and  cause  a  disturbance-dependent  correcticxi  to  the  stability  limit  defined  by  small  per¬ 
turbation  techniques. 

Evidence  of  the  occurrence  of  this  phenomenon  in  the  present  case  is  provided  by  the 
contrast  between  Figures  23.1  b,  b,  c  and  figures  39  a,  b.  Both  sets  of  figures  refer  to  Run 
6  with  B  =  10  and  the  same  initial  position  and  velocity.  However,  in  the  Figure  23.1  series, 
the  initial  pressure  distribution  is  closer  to  the  steady  state  than  in  the  Figure  39  series. 

As  the  figures  show  and  the  numerical  results  confirm,  the  orbit  of  Figures  24.1  a,  b,  c  (smal¬ 
ler  perturbation)  is  convergent  and  the  orbit  of  Figures  39  a,  b  (larger  perturbation)  is 
divergent. 


The  above-mentioned  reasons  contributed  to  the  formati<»  of  the  following  policy  in 
the  running  of  finite  orbit  cases;  the  threshold  value  of  B  was  defined  within  limits  which 
are  narrow  enough  for  practical  purposes  but  not  further;  small  perturbations  of  ecjuilibrium 
c(xiditions  were  used  in  order  to  provide  a  check  on  the  accuracy  of  the  linearized  theory  of 
Section  IV.  A  study  of  the  non-linearity  effects  aimed  at  developing  sufficioit  stability 
criteria  on  a  quantitative,  rather  thn  c|ualitative  basis  was  not  undertaken,  because  presently 
outside  the  mesns  of  the  project. 

From  a  c^alitative  point  of  view  it  can  be  aaid  that  atable  operation  waa  alwaya  ob¬ 
served  for  values  of  B  larger  than  a  threahold  10  to  15%  higher  than  that  predicted  by  the 
linearized  theory.  In  using  such  numbers,  however,  cne  should  reaMuber  that  no  systematic 
study  of  this  effect  has  been  made,  and  that  these  limits  might  well  be  exceeded  for  some 
combination  of  ruining  parameters. 
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One  more  interesting  non-linear  effect  has  not  been  studied  in  detail;  the  occur¬ 
rence  of  limit  cycle  or  "  stable  orbits".  Experimental  evidence  indicates  that  stable  finite 
orbits  actually  occur  in  rotating  machinery.  For  long  complete  gas  journal  bearings  evidence 
is  not  plentiful  because,  after  the  inception  of  half- frequency  vdiirl,  expensive  bearing  fail¬ 
ure  usually  occurs  almost  immediately.  Actual  machines  are  never  completely  rid  of  unbalance 
and  bearing-surface  out-of-roundneas,  so  that  dynamic  stabilizing  effects  from  these  agencies 
(in  the  form  of  squeeze  film  forces  snd  phase  angle  shifts)  can  help  establish  an  orbital 
equilibrium  which  would  be  impossible  for  a  perfect  shaft.  (X>  the  other  hand,  it  can  be  very 
expensive  to  make  a  systematic  study  of  limit  cycles  by  means  of  the  orbit  program  because  it 
is  very  difficult  to  distinguish  between  the  tendency  toward  an  orbit  frcm  simple  slow  con¬ 
vergence  or  divergence. 

Ihe  results  of  the  above-mentioned  runs  designed  to  check  the  accuracy  of  the  small 
perturbation  progrmn  are  condensed  in  Table  III.  For  the  purpose  of  better  visualization  these 
results  have  been  plotted  on  stability  maps  of  the  type  of  Figure  4.  For  the  sake  of  clarity 
the  stability  map  has  been  split  in  two;  Figure  23  for  6^  =  0.2,  0.6,  0.9  and  Figure  24  for 
6^  =  0.4,  0.8.  Ihe  points  corresponding  to  a  given  case  lie  on  a  A  =  constant  line  and 
are  represented  by  a  triangle  if  the  orbit  is  stable  and  a  dot  if  the  orbit  is  unstable. 

A  series  of  orbits  as  obtained  from  the  CHARACTRCN  is  shown  in  Figures  25  through 
34.  Unfortunately,  the  scaling  routine  uses  a  minimum  plot  size  of  0.2  x  0.2  in  X  and  Y  so 
that  the  orbits  representing  the  system  response  to  small  disturbances  appear  as  smeared  dots 
and  conclusions  on  stability  have  to  be  drawn  from  the  numerical  list  of  X-Y  pairs.  Such  a 
situation  is  shown  by  Figures  27  a,  b,  and  c  and  is  true  for  all  points  which  are  close  to 
the  stability  threshold.  The  reader  will  remember  that,  within  certain  limits,  the  disturbance 
size  plays  a  role  in  the  definition  of  instability  and  therefore  will  understand  why  such 
small  orbits  were  necessary.  Since  no  definition  is  offered  by  such  plots,  only  one  sample 
run  is  shown.  Figures  25  a  through  41  b  are  illustrations  of  orbits  with  the  particular  fea¬ 
tures  discussed  in  this  section. 


Figures  25a  through  34  b  are  illustrations  of  stable  and  unstable  orbits  for  Runs 
7,  8  and  12. 

Fig  ires  35  a  through  37  b  are  orbits  corresponding  to  Runs  17  and  16  and  show  the 
long  duration  of  transients  in  cases  of  high  values  of  A. 


Figures  38  a  and  b  show  the  response  of  Rin  17  when  an  arbitrary  amount  of  unbalance 
mass  Uj  is  attached  at  a  radius  Rj  from  the  shaft  center  so  that 

4^1  R.  =  , 


Me 


The  striking  feature  of  this  orbit  is  that  it  is  very  similar  to  one  which  was  experimentally 
obtained  by  Elwell  [12]  of  the  General  Electric  (^pany.  This  orbit  should  be  interpreted 
only  from  a  qualitative  standpoint  since  no  systematic  study  of  the  effect  of  unbalance  has 
been  carried  out. 


Figures  28  a,  b,  39  a,  b,  and  40  b,  are  illustrations  of  possible  limit  cycles  aris¬ 
ing  from  two  of  the  attached  runs. 

Figures  41  a,  b,  show  an  orbit  which  has  grown  to  cover  nearly  the  entire  clearance 
area.  It  can  be  seen  that  the  divergence  rate  decreases  once  the  amplitude  of  oscillation  be¬ 
comes  large  enough  to  boost  the  effect  of  squeeze  film  forces. 
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V.  DISCUSSION 


1.  General  Diacu— ion  md  Ceupwiacn  with  the  of  Oth«r  Inve«tit*tw 

‘Die  eucceee  of  the  t«N)  techniquee  presented  in  this  paper,  aid  the  results  they 
have  produced,  vdll  help  settle  a  nimher  of  questions  «4iich  have  been  pussling  research  engi¬ 
neers  for  some  tiise. 

It  was  previously  mentioned  that  one  of  the  nqjor  obstacles  encountered  in  the  study 
of  gas  bearing  dyntmic  problems  has  been  the  treatment  of  the  "history”  effect.  Obviously, 
this  effect  prevents  the  use  of  field  charts  of  gas  film  forces  expressed  as  instantaneous 
functions  of  position  and  velocity.  IhMiks  to  the  orbit  pressures,  it  is  now  possible  to 
evaluate  the  relative  ma^iitude  of  terms  in  Reynolds'  eqiation.  Figures  42a  throu|^  42k,  and 
43,  and  44  all  concern  Rin  17  with  B  =  10  and  contain  respectively:  the  ffl  distributions  at 
every  time  step,  the  values  of  X,  Y,  t,  at  every  tisie  step,  and  the  orbit  covered  by  300 
steps.  It  will  be  of  interest  to  evaluate 
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Some  sample  calculations  for  the  point  6*0*  (last  PH  point)  and  between  steps  6103 
and  6118  yield  the  results  listed  in  Table  IV.  It  crni  be  easily  seen  that  the  orders  of  magni¬ 
tude  of  H  ^  and  P  ^  are  conqiarable,  and  no  strict  justification  exists  for  neglecting 

either  term.  More  pressure  distributions  are  available  both  from  this  run  aid  from  other  ones 
so  that  more  extensive  studies  can  be  conducted  to  evaluate  the  ma^iitude  of  the  serious  tenas 
in  Reynolds’  equation  in  various  orbital  conditions. 


It  should  be  pointed  out  that  this  case  was  chosen  at  random  and  that  the  relative 
importance  of  these  terms  might  vary  in  other  cases. 
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An  iaportant  fact  to  be  noticed  in  all  orbita  obtained  up  to  the  preaent  tine  ia  that 
the  journal  canter  path  aeana  to  have  little  tendency  to  encloae  the  bearing  center.  Ihie  in- 
fomation  night  throe  unfavorable  light  on  approxinetione  ehich  aaaune  etea^  orbiting  of  the 
journal  canter  at  alnoat  conatnt  eccentricity  nd  ingular  velocity. 

Conparisona  with  other  theoriee  ahow  that  beat  agreenent,  both  in  the  qualitative 
and  quantitative  aenae,  ia  achieved  by  Qieng'a  analyaia  [S].  Ihia  waa  expected,  becauae  Galer- 
kin'a  method  yielda  a  aurpriaingly  good  ^>proximation  to  ateady-atate  preasure  diatribution. 

Ihe  quaai-atatic  thewy  of  Stemlicht  and  Rentaepia  [26]  cornea  relatively  close  to  the  anall 
perturbation  reaulta  but  it  ia  quite  inadequate  at  low  valuea  of  C*  (high  values  of  A).  Ibis 
may  be  understood  when  we  realize  that  neglect  of  compressibility  effects  will  lead  to  greater 
errors  at  high  valuea  of  the  ccmpressibility  parameter  A. 

Auaman's  reaulta  ahow  great  discrepancy  with  all  other  theories.  This  fact  can  be 
attributed  to  two  basic  reasons:  first,  hia  linearized  ffl  expansion  does  not  give  the  high  ac¬ 
curacy  needed  in  stability  calculations  even  at  moderate  values  of  eccentricity  ratio  sec¬ 
ond,  for  the  purpose  of  mathematical  expediency,  the  terms 


1  - 


X2  ♦  >1 1  -  Y 


1  - 


>1 1  - 


and  ( 1  -  vl  1  -  ^ 


are  arbitrarily  taken  to  be  constant  in  time.  It  is  also  possible  that  some  errors  in  calcula¬ 
tions  have  been  made  in  the  evaluation  of  the  preliminary  stability  maps  received  by  the 
authors. 


2.  Practical  Value  of  this  Investigation  and  Sample  Problem 

It  is  well  recognized  in  practice  and  it  has  been  the  result  of  all  theories  that, 
for  the  same  dimensionless  load,  the  shorter  the  bearing  the  higher  the  instability  threshold 
speed.  This  effect  is  amenable  to  intuitive  justification  when  we  consider  that,  for  a  short 
bearing,  the  same  load  produces  larger  eccentricities  ratios  in  steady  state  operation.  Then 
the  squeeze-film  effects  produce  both  higher  overriding  response  frequencies  md  larger  danp- 
ing  due  to  axial  inflow  and  outflow.  By  the  reasoning  presented  in  Section  V.  4,  these  effects 
correspond  to  higher  threshold  speeds.  Therefore,  from  this  practical  point  of  view,  the 
results  of  any  accurate  theory  based  on  bearings  of  infinite  axial  length  give  conservative 
stability  criteria. 

Unfortunately,  a  pessimistic  criterion  in  the  present  stability  problem  is  often  not 
adequate  for  practical  applications.  The  tendency  to  whirl  is  so  pronounced  that  using  pessi¬ 
mistic  criteria  could,  in  many  applications,  rule  out  a  considerable  portion  of  the  useful 
range  of  operating  parameters.  However,  extension  of  the  treatments  presented  in  this  paper 
to  finite-length  bearings  is  not  expected  to  present  any  difficulties  in  principle. 

The  orbit  program  can  be  very  useful  for  practical  applications  since  it  takes  into 
account  non-linearity  and  will  give  valid  answers  for  real  situations.  Extension  of  this  pro¬ 
gram  to  the  finite- bearing  case  is  being  carried  out  at  Columbia  University  by  the  author. 

A  typical  design  problem  for  a  long  bearing  will  require  the  determination  of  the 
clearance  ratio  (VR  which  will  allow  stable  operation  of  a  bearing  carrying  a  dimensionless 
load  L  at  an  eccentricity  ratio  not  larger  than  £  and  at  rotational  speed  H. 
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The  conbination  of  minia»  lood  and  ■Mcina  e  lead  to  the  apecification  of  a  aini* 
nun  A  from  Elrod  and  Burgdorfer  reaulta  of  Figure  2.  Having  evaluated  (tf  fron  Imoan  phyaical 
paranetera  ind  Q,  «e  can  enter  Fig«re  5  and  deteradne  the  m vinn  C*  correapondiag  to  A  and 
of.  The  converaion  from  the  naxintm  value  of  C*  to  the  nexinaim  value  of  C/R  ia  obvioua  from 
the  definition  of  C. 

Aa  a  nuii«rical  exanple  let  ua  conaider  an  air  bearing  to  operate  with  a  unit  load 
L  =  2  at  e  =  0.8 
From  Figure  2  this  corresponds  to 

A  =  0.7 

Then 


-  4.2  from  Figure  4 


or 


B 


=  0.113 


We  can  new  compute 

,0.2 


C*  = 


(^) 


=  2.05 


If  the  load  was  due  to  gravity 


\7^2  I  R  \72  g  ” 


Then,  for  R  -  2  in, 
C  =  1.336  X  10^ 


which  yields 


(£)<  i: 

'^'stable 


05 

336 


10 


•3  - 


From 

6  jj.  CO  V? 

A  =  -  ,  CO 

P.C? 


1.53  X  10*^ 


AP.C„2 


U90  —  ~  14,200  mi 
sec 
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3.  Gancluding  Ranarks  id  Bec<ininend»tions 

Ihe  experience  of  ell  experimental  investigators  in  the  gas  bearing  field  has  been 
that  inat^ility  ia  the  major  problem  to  overccme.  Stability  maps  for  the  idealized  model 
treated  in  this  paper,  even  with  eventual  extension  to  finite- length  bearings,  certainly  repre¬ 
sent  a  forward  step  in  understanding  the  mechaniam  of  this  phenomenon,  but  should  not  be  inter¬ 
preted  as  the  final  answer.  In  fact,  strictly  speaking,  the  dynamics  of  the  shaft  in  its  bear¬ 
ings  cannot  be  studied  separately  from  the  rest  of  the  system.  Bearing  supports,  machine  frame, 
driving  mechanisms,  unbalance,  shaft  flexibility  are  elements  which  operate  in  direct  coupling 
with  the  rotor-bearing  system.  Obviously,  it  would  be  extremely  costly  and  so^owhat  wasteful 
to  develop  a  general  theory  taking  all  factors  into  account.  It  is  possible,  however,  to  con¬ 
struct  computer  programs  similar  to  the  orbit  progran,  but  generalized  to  include  the  above- 
mentioned  effects.  Runs  would  then  be  made  only  for  the  purpose  of  aiding  the  design  of  parti¬ 
cular  machines.  Ihis  tool  would  be  of  extreme  practical  value  because  their  low  level  of  internal 
damping  makes  gas  bearings  susceptible  to  harmful  resonances  with  natural  frequencies  of  some 
other  members  of  the  machine  structure. 


With  the  present  stability  theories  now  at  hand,  it  would  seem  quite  appropriate  that 
extensive  sets  of  accurate  experiments  be  run  to  give  a  firmer  footing  for  further  theoretical 
efforts.  Furthermore,  instrumentation  and  machining  techniques  have  now  been  developed  ahich 
are  capable  of  handling  physical  dimensions  in  microinches  with  a  fair  degree  of  reliability, 
so  that  such  experimentation  is  definitely  possible. 

<&dtLr 
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NOMENCLATURE 


Dbfinition 


R  Pj  (tC  ~~)  Dunentionleas  Stability  Parameter 
4 

Radial  Qearance 

(2/B^2)l/5  =  •  Q/R  -  Dimenaionleaa 

critical  clearance  Paraneter 


eccentricity 

Pressure  Force  Ccmponents  in  x,  y  direction 

Dimensionless  Film  Ihickness 

Film  Thickness 

Function  of  s 

Function  of  s 

Function  of  s 

Friction  Force  Components  in  x,  y  direction 
W/(R  .  P^)  Diaiensionless  Load  Parameter 
Function  of  s 

Rotor  Mass  per  Unit  Axial  Length 
Film  Pressure 
Anbient  Pressure 


—  Dimensionless  Pressure 

P. 


Function  of  0  and  s 


AT 


Nunerical  Stability  Parameter 


Shaft  Radius 
Function  of  0  and  S 
Viscous  Shear  Stresa 


(%.  of  first 
apfiearance 

8 

4 

20 

6 

4 

4 

18 

18 

18 

7 

8 

18 

8 

4 

4 

4 

12 

24 

4 

12 

7 
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T=  ((V2)‘ 

t 

W 

X  =  Xj/c 

Y  =  Yj/c 

€ 

T)  =  /3/R 

V 

e 

A 

M 


PH 

n 

A 

CO 

Cl* 

of 

(  )“ 

(  )' 

(  -  ) 


NOMENCLATURE  (Coat.) 

Definition 

Tmefonaed  Time  Variable 
Dimenaionless  Time  Parameter 
Time 

External  Load  per  IMit  Axial  Length 
Dimensionless  Cartesimi  Qx>rdinate  of  Shaft  Center 
Dimensionless  Cartesian  Coordinate  of  Shaft  Center 
Coordinate  Nomal  to  Direction  of  Motion 
Eccentricity  Ratio  =  e/c 

Dimenaionless  Luigth  Normal  to  Direction  of  Motion 
Dimensionless  Coordinates 

Angle  Measured  in  the  Direction  of  Journal  Rotation 
(£1  ^  R^/(P^C^)  Dimensionless  Speed  Parameter 
Dynamic  Viscosity  Coeff. 

K  'll  -  c* 

^  ~A 

Dimensionless  Coordinate 
Dimensionless  Variable 
Angular  Velocity 

MCr^/(^RP^)  =  2/('rrB)  =  Dimensionless  Stability  Parameter 

Complex  part  of  s 

>1 4/{LB)  =  Stability  Parameter 
3  u  M*  1/5 

(WB^)  1/5  =  (“ — ry)  Dimenaionless  Critical 

^  Speed  Parameter 

Parameter  or  Coordinate  at  Equilibrium  Condition 

Differentiation  with  Respect  to  6 

LaPlace  Transform  of  a  Fuiction 
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Pg.  of  first 
appearance 

11 

4 

4 

4 

4 

4 
7 

5 
5 

4 

5 

4 

19 

22 

5 
4 

20 
12 
20 
20 

9 

10 
11 
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belli 


Am 

V 

A 

L=  2C^ 

Yo 

1 

1.0 

.2 

.9622 

.4414 

.14326 

.13956 

2 

.4 

.8209 

.8312 

.29710 

.26782 

3 

.6 

.5418 

.9036 

.47781 

.36290 

4 

.8 

.2052 

.6078 

.65156 

.46418 

5 

.9 

.0722 

.3416 

.71306 

.54914 

6 

0.5 

.2 

1.979 

.5828 

.09059 

.17831 

7 

.4 

1.859 

1.2734 

.18840 

.35285 

8 

.6 

1.460 

1.9934 

.30662 

.51574 

9 

.8 

0.6348 

1.9074 

.45943 

.65492 

10 

.9 

0.2221 

1.1866 

.48334 

.75920 

11 

0.25 

.2 

4.0126 

0.6428 

.04850 

.19403 

12 

.4 

3.9754 

1.5068 

.09936 

.38746 

13 

.6 

3.643 

2.880 

. 15630 

.57928 

14 

.8 

2.379 

4.894 

.23701 

.76409 

15 

.9 

1.017 

4.596 

.30487 

.84679 

16 

.125 

.2 

8.058 

.6604 

.02484 

.19845 

17 

.4 

8.110 

1.5738 

.04990 

.39688 

18 

.6 

7.781 

3.162 

.07641 

.59512 

19 

.8 

6.256 

6.758 

.10673 

.79285 

20 

,9 

4.246 

10.922 

.12940 

.89064 

21 

0.0 

1.0 

0.0 

0.0 

0.0 

22 

1.0 

.1 

.99107 

.221950 

.07094 

.07048 

23 

.5 

.1 

1.9957 

.283682 

.04485 

.08937 

24 

.25 

.1 

4.0043 

.30896 

.02428 

.09701 

25 

.125 

.1 

8.01657 

.31633 

.01242 

.09923 
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TABLE  1  A 

INSTABILITY  THBEaHOLDS 


An 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
22 

23 

24 

25 


B 

9.4 

2.45 

0.97 

0.36 

0.116 

5.80 

1.2 

0.47 

0.115 

0.0021 

5.6 

1.15 

0.37 

0.114 

0.0045 

5.0 

1.1 

0.345 

0.11 

0.038 

32. 

21. 

19. 

18. 


"1* 

0.96 

1.40 

2.14 

4.27 

10.05 

1.09 

1.62 

2.07 

4.27 

40.1 

1.06 

1.52 

1.94 

2.68 

13.9 

1.10 

1.52 

1.91 

2.32 

3.11 

.75 

.819 

.826 

.838 
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TABLE  II 

COMPARISON  OF  EQUILIBRIiai  DISTOIBUTIOWS* 


Elrod  &  BUrgdorfer  Caatelli  &  Elrod  Oirbit  ProgrM 


d  (deg.) 

PH 

d  (deg.) 

PH 

194.1 

0.93435 

12 

1.30588 

187.6 

0.96949 

24 

L 39255 

180.6 

1.00953 

36 

1.46353 

173.1 

1.05558 

48 

1.51551 

164.8 

1. 10892 

60 

1.54646 

155. 5 

1.17090 

72 

1.55562 

144.9 

1.24248 

84 

1.54344 

132.8 

1.32331 

96 

1.51145 

118.6 

1.40969 

108 

1.46217 

102.0 

1.491% 

120 

1.39889 

82.8 

1.54807 

132 

1.32547 

61.0 

1.55058 

144 

1.24602 

37.5 

1.47433 

156 

1.16463 

13.7 

1.32218 

168 

1.08502 

350.8 

1. 13351 

180 

1.01030 

330.1 

0.96086 

192 

0.94271 

312.0 

0.83789 

204 

0.88364 

296.5 

0.76765 

216 

0.83365 

283.2 

0.73570 

228 

0.79279 

271.8 

0.72637 

240 

0.76088 

261.8 

0.72927 

252 

0.73804 

253.0 

0.73884 

264 

0.72523 

245.1 

0.75223 

276 

0.72479 

237.8 

0.76811 

288 

0.74034 

231.0 

0.78593 

300 

0.77572 

224.6 

0.80544 

312 

0.83288 

218.5 

0.82669 

324 

0.91034 

212.4 

0.64983 

336 

1.00331 

206.4 

0.87517 

348 

1.10498 

200.3 

0.90315 

360 

1.20810 

*Caae  of 


V  =0.5 

6o  =  0.6 

A  =  1.46 
L  =  1.9934 
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TABLE  III 


Rm  AT  y. 

o  o 


"I* 


Stable? 


7  .026  at. at. 


8  .011  at.  at. 


12  .026  at. at. 


13  ,  026  at.  at. 


16  .026  St. at. 


17  0.026  St. St. 


18  .026  at. at. 


19  .026  St. at. 


at. St.  0  0 


St.  at.  0  0 


at. St.  0  0 


St. St.  0  0 

.01  .01 

.01  .01 

0  0 

0  0 

St. St,  0  0 


St. St.  0  0 


St. St.  0  0 


St. St.  0  0 


8 

.627 

yes 

4 

.886 

yea 

1.8 

1.321 

yea 

1.3 

1.555 

yes 

1.1 

1.690 

no 

.8 

1.982 

no 

4 

.m 

yea 

.8 

1.584 

yes 

.5 

2.004 

yea 

.45 

2.112 

no 

.2 

3.167 

no 

1.5 

1.330 

yes 

1.2 

1.487 

no 

1.0 

1.629 

no 

7.0 

.445 

yes 

0.6 

1.521 

yes 

0.45 

1.757 

yea 

0.3 

2.152 

no 

0.07 

4.454 

no 

20 

.550 

yes 

8 

.870 

yes 

7 

.930 

no 

5 

1.101 

no 

1 

2.461 

10 

.504 

yea 

7 

.603 

yea 

6 

.651 

yea 

1 

1.59 

no 

0.3 

2.91 

no 

8 

.398 

yea 

4 

.562 

yes 

.8 

1.258 

yes 

.35 

1.901 

no 

.25 

2.250 

no 

0.14 

1.056 

yea 

0.11 

1.320 

no 

0.08 

2.721 

no 

0.07 

2.288 

yen 

.04 

3.026 

no 

.02 

4.280 

no 

20 


026 


at.  at. 


St.  St. 
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TABLE  IV 


Step 

(**  n)/\  it) 

6103 

.1.5635 

6104 

-1.5187 

6105 

-1.5012 

6106 

-1.5069 

6107 

-1.5377 

6108 

-1.6022 

6109 

-1.7278 

6110 

-1.9836 

6111 

-2.6635 

6112 

-8.0809 

6113 

♦1.9262 

6114 

♦0.2063 

6115 

-0.2307 

6116 

•0.4238 

6117 

-0.5282 

6118 
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CASE  NUMBER  \i  e>  O.UttOOO' 


STEP 

X 

Y 

STEP 

X 

Y 

1 

0.1574864 

0.5795248 

51 

0.  1  (.28fci35 

0.5813254 

2 

0.1577459 

0.5797841 

52 

O.li.2  7789 

0.5810893 

3 

0.1580045 

0.5800420 

53 

0.1 (26676 

0.5808479 

4 

0.1582617 

0.5802976 

54 

0.1' 25498 

0.5806020 

S 

0.15851 70 

0.5805502 

55 

0. 1124258 

0.5803520 

6 

0. 1587700 

0.5807991 

56 

0. 1 122959 

0.5800987 

7 

0.1590202 

0.5810436 

57 

0.112 1605 

0.5798428 

8 

0.1592672 

0.5812830 

58 

0.1 120198 

0.5795848 

9 

0.159S105 

0.5815166 

59 

0. M 18741 

0.5793254 

10 

0.1 597497 

0.5817437 

60 

0.  1..  1  7238 

0.5790654 

n 

0.1599845 

0.5819637 

61 

0. 1ii15693 

0.5788053 

12 

0.1602143 

0.5821 760 

62 

0.  Ill  14108 

0.5785456 

13 

0. 1604389 

0.5823800 

63 

0.  11112486 

0.5782876 

14 

0.1606578 

0.5825750 

64 

0.  I(il0832 

0.5780313 

IS 

0.1608707 

0.5827606 

65 

0.1609148 

0.5777775 

16 

0.1610772 

0.5829363 

66 

0. 1607438 

0.5775269 

17 

0.1612770 

0.5831014 

67 

0.1635705 

0.5772600 

18 

0.1614698 

0.5832557 

68 

0. 1603952 

0.5770374 

19 

0.1616552 

0.5833966 

69 

0.1602184 

0.5767999 

20 

0.1618331 

0.5835298 

70 

0. 1600402 

0.5765676 

21 

0.1620030 

0.5836490 

71 

0. 1578610 

0.5763417 

22 

0.1621649 

0.5837557 

72 

0.1576812 

0.5761223 

23 

0.1623184 

0.5838497 

75 

0.1575011 

0.5759099 

24 

0.1624633 

0.5839308 

74 

0. 1593209 

0.5757052 

2S 

0.1625995 

0.5839987 

75 

0.1591410 

0.5755084 

26 

0.1627267 

0.5840533 

76 

0.1589616 

0.5753202 

27 

0.1628448 

0.5840944 

77 

0.158/831 

0.5751408 

28 

0.1629537 

0.5841220 

78 

0.1586058 

0.5749708 

29 

0.1630532 

0.5841359 

79 

0.1584298 

0.5746*04 

30 

0.1631433 

0.5841362 

80 

C. 1582556 

0.5746601 

31 

0.1632239 

0.5841228 

81 

0.1580833 

0.5745200 

32 

0.1632949 

0.5840959 

82 

C. 15/9131 

0.5743906 

33 

0.1633563 

0.5840554 

83 

C. 15/7454 

0. 5742720 

34 

0. 1634080 

0.5840016 

84 

C. 15/5803 

0.5741645 

35 

0. 1634501 

0.5839346 

85 

C. 1574 181 

0.5740683 

36 

0.1634826 

0.5838545 

86 

0.1572590 

0.5739836 

37 

0.1635055 

0.5837617 

87 

0.1571031 

0.5739105 

38 

0. 1635189 

0.5836564 

88 

0.1569507 

0.5738490 

39 

0.1635228 

0.5835389 

89 

0.1568019 

0.5737994 

40 

0. 1635174 

0.5834095 

90 

0.1566569 

0.5737616 

41 

0.1635028 

0.5832685 

91 

0 . 1 565158 

0.5737356 

Ffc,. 

42 

0.1634790 

0.5831165 

92 

0.1563789 

0.5737215 

43 

0.1634462 

0.5829538 

93 

0.1562461 

0.5737193 

13 

44 

0.1634047 

0.5827808 

94 

0. 1501176 

0.5737287 

45 

0.1633545 

0.5825980 

95 

0.1559935 

0.5737498 

ir  mo.ozg 

46 

0. 1632958 

0.5824059 

96 

0.1558740 

0.5737824 

47 

0.1632290 

0.5822051 

97 

0.1557590 

0.5736264 

Olf'ftlT 

48 

0.1631540 

0.5819961 

98 

0.155648/ 

0.5736816 

49 

0.1630713 

0.5817794 

99 

0.1555431 

0.5739477 

50 

0.1629811 

0.5815557 

100 

0.1554423 

0.5740246 
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CASE  NUMBER  1i  B«  O.USOOO 


step 

1  0.1 

2  0.1 

\  O.I 
4  0.1 

b  0.1 
6  O.I 
'  0.1 
U  0.1 
9  0.1 

10  0.1 

11  O.I 

12  0.1 
li  0.1 

14  0.1 

15  O.I 

16  O.I 

1  /  0.1 

18  O.I 

19  0.1 

20  O.I 

21  0.1 

22  O.I 

23  0.1 

24  0.1 

25  0.1 

26  O.I 
21  0.1 
28  0.1 

29  0.1 

30  O.I 

31  O.I 

32  O.I 

33  0.1 

34  0.1 

35  O.I 

36  O.I 
if  0.1 

38  O.I 

39  O.I 

40  O.I 

41  O.I 

42  O.I 

43  O.I 

44  0.1 

45  O.I 

40  0.1 

47  O.I 

48  O.I 

49  O.I 

50  O.I 


X 

573564  0 

574863  -  0 
576162  0 

57/458  0 

578752  0 

580043  C 
58  1331  0 

5o2ol5  0 
583894  0 

585167  0 

586435  0 

587697  0 

S88951  0 

590198  0 

591437  0 

59266/  0 

593889  0 

595100  0 

590 302  0 

59/492  0 

598672  0 

599840  0 

600995  0 

602138  0 

603268  0 

604384  0 

605486  0 

6065  73'  0 

607646  0 

608702  0 

609743  0 

610768  0 

611776  0 

612767  0 

613740  0 

614695  0 

615632  0 

616551  0 

617450  0 

618330  0 

619191  0 

620031  0 

620851  0 

621651  0 

622430  0 

623187  0 

623924  0 

624638  0 

625331  0 

626002  0 


Y 

.5/93948 
.5795247 
.5796544 
.5797839 
.5799130 
.5800416 
.5801697 
.5802971 
.5804238 
.5805497 
.5806746 
.5807985 
.5809213 
.58  10429 
.5811633 
.5812822 
.5813998 
.5815158 
.5816302 
.5817429 
.5818538 
.5819628 
.5820700 
.5821751 
.5822781 
.5823790 
.58247/7 
.5825741 
.5826681 
.5827597 
.5828488 
.5829355 
.5830193 
.5831006 
.5831791 
.5832549 
.5833278 
.5833979 
.5834650 
.5835291 
.5835903 
.5836483 
.5837033 
.5837551 
.5838038 
.5838492 
.5838915 
.5839304 
.5839661 
.5839985 


STEP 

51 

52 

53 

54 

1!  . 

55 

,) 

56 

0. 

57 

.1. 

58 

(.  . 

59 

.1. 

60 

0. 

61 

62 

0. 

63 

0. 

64 

0  . 

65 

0. 

66 

0. 

67 

0. 

68 

0. 

69 

0. 

70 

0. 

71 

0. 

7? 

0 . 

73 

0. 

74 

0. 

75 

0. 

76 

0. 

77 

5. 

78 

l) . 

79 

3. 

bO 

). 

bl 

3. 

82 

■3  . 

83 

0. 

84 

0. 

H5 

^  ♦ 

86 

0. 

o7 

3  . 

88 

3. 

89 

0  • 

90 

u  • 

91 

3. 

92 

0 . 

93 

0. 

94 

0. 

95 

3. 

96 

3. 

97 

3. 

98 

A 

M  • 

99 

0. 

100 

3. 

X 

1  .5;i6650 
15;;  72  76 
15;-7879 
1 5;;3459 
15;-90  17 
1  i;-955i 
I  .>::336  I 
1530549 

15  110  12 
16.11452- 
16:11868 
16:1226  I 
16.12829 

16  12974 
16.15294 
16.13590 
1613863 
16.141  1  I 
1634335 
16.14535 
16.1471  I 
16.14864 
16.14992 
16.15096 
16  151  77 
1615234 
1615267 
16 152/7 
16.15263 
16*5226 
16.15166 
l6  <5084 
16149/8 
16 14850 
16.14699 
16 14526 
16  14.331 
16*4  1  15 
16.138  76 
1633616 
1613336 
16*3034 
16  12712 
1612369 
1632006 
1631624 
1611222 
1630801 
1630360 
1629902 


Y 

0.5840275 
0.5840532 
0.5840755 
0.5840944 
0.5841 100 
0.584122  I 
0.5841 309 
0.5841362 
0.5841381 
0.5841 366 
0.5841 317 
0.5841234 
0.5841117 
0.5840966 
0.5840781 
0.5840562 
0.584031 1 
0.5840025 
0.5839707 
0.5839356 
0.5838973 
0.5838557 
0.5838109 
0.5837630 
0.5837119 
0.583657b 
0.5836005 
0.5835403 
0.5834771 
0.5834  109 
0.5833419 
0.5832701 
0.5831954 
0.5831  ISO 
0.5830380 
0.5829553 
0.5828701 
0.5827823 
0.5826921 
0.5825995 
0.5825046 
0.5824074 
0.5823080 
0.5822065 
0.5821029 
0.5819973 
0.5818898 
0.581 7805 
0.5816694 
0.58155^6 


Fic^.  2a  cu 
/3 

0.01^ 


ORSiT 

C  OO  CD'  W  ATE  i 
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CASE  NUMBER  It  R*  0.4S00I 


TEP 

X 

y 

STEP 

;< 

Y 

101 

0.162942S 

0.S8 14421 

151 

C  -  15911621 

0.5754002 

102 

0. 1628930 

0.S8I326I 

15? 

C. 15M"723 

0.5753068 

lOi 

0.1628417 

0.S8l20a7 

153 

C. 1680828 

0.5752156 

104 

0.1627888 

0.S810898 

154 

C. 1587634 

0.5751267 

10S 

0.1627341 

0.S809696 

155 

C. I587C44 

0.5750401 

106 

0.1626777 

0.S808482 

156 

C.  l5ili>IS6 

0.5749559 

107 

0.1626198 

0.S8072S6 

157 

0. I58S272 

0.5748  741 

108 

0.162S603 

0.S8060i9 

158 

0.  15(li-39t 

0.5747947 

109 

0.1624992 

0.S804772 

159 

0.1585515 

0.5747179 

110 

0.1624366 

0.S803S16 

160 

0.1582643 

0.5746436 

111 

0.162372S 

0.S8022S2 

161 

0.1581776 

0.5745719 

112 

0.1623070 

0.S800980 

162 

0.1580914 

0.5745028 

113 

0.1622401 

J.S79970I 

163 

0. 158C057 

0.5744363 

1 14 

0.1621718 

0.S7984I6 

164 

0. 157??06 

0.5743726 

ns 

0.1621022 

0.S797I27 

165 

0.157&36I 

0.5743116 

1 16 

0.1620314 

0.S79S832 

166 

0.1577522 

0.5742533 

117 

0.I619S93 

0.S794S3S 

167 

0.1576690 

0.5741978 

118 

0.I6I88S9 

0.S793234 

168 

0.1575864 

0.5741451 

119 

0.16I8I1S 

0.S791932 

169 

0.1575045 

0.5740953 

120 

0.161 73S9 

0.S790629 

170 

0. 1574234 

0.5740482 

121 

0.1616S92 

0.S789326 

171 

0.1573430 

0.5740041 

122 

o.i6isais 

0.S788023 

17? 

0.1572635 

0.5739629 

123 

0.I61S028 

0.S786722 

173 

0.1571847 

0.5739245 

124 

0.1614231 

0.S78S423 

174 

0.1571067 

0.5738891 

I2S 

0.1613426 

0.S784127 

175 

0.1570296 

0.5738566 

126 

0.1612611 

0.S78283S 

176 

0.1569534 

0.573R271 

127 

0.1611788 

0.S781S48 

177 

0.1568781 

0.5738005 

128 

0.16109S8 

0.S780266 

178 

0.1568036 

0.5737769 

129 

0.1610120 

0.S778991 

179 

0.1567302 

0.5737563 

130 

0.1609274 

0.S777722 

180 

0.1566577 

0.5737386 

131 

0.1608423 

0.S776461 

181 

0 .1565861 

0.5737239 

132 

0.I607S6S 

0.577S209 

18? 

0.1565156 

0.5737122 

133 

0.1606701 

0.S7 73966 

183 

C. 1564460 

0.5737035 

134 

0.160S83I 

0.S772734 

184 

0.1565775 

0.5736977 

I3S 

0.I6049S7 

0.S771S12 

185 

C. 1565100 

0.5736949 

136 

0.1604078 

0.5770301 

186 

0. 1562436 

0.5736950 

137 

0.I60319S 

0.S769I03 

187 

0. 1561783 

0.5736982 

138 

0.1602309 

0.S7679IH 

186 

0. 1561 140 

0.5737042 

139 

0.160  1419 

0.S766747 

189 

0.1560508 

0.5737132 

140 

0.1600S26 

0.5765590 

190 

0. 1559388 

0.5737251 

141 

0.1S99630 

0.5764449 

191 

0.  1559278 

0.5737398 

142 

0.1698732 

0.5763323 

19? 

0. I55do80 

0.5737575 

143 

0.1S97833 

0.5762213 

193 

0.  1558.194 

0.5737780 

144 

0.IS96932 

0.576II2I 

194 

0. 1557518 

0.5738013 

14S 

0.IS96030 

0.5760046 

195 

0.1556955 

0.5738275 

146 

0.1S9SI28 

0.5758989 

196 

0. 1556403 

0.5738565 

147 

0.1S9422S 

0.5757952 

197 

0.1555863 

0.5738882 

148 

0.IS93325 

0.5756934 

198 

0. 1555334 

0.5739226 

149 

0.1S92421 

0.5755936 

199 

0.1554818 

0.5739598 

SO 

0.lS9iS2l 

0.5754958 

200 

0. 1554313 

0.5739996 

rio,.  20  h 

*  0.013 


O  <B>T 
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APPENDIX  1 

DERIVATION  OF  LUBRICATION  EQUATIONS 


Let  u,  V,  w  be  the  velocity  V  component  in  the  a,  /S,  y  direction!  respectively.  Then 
the  equation  of  continuity  is 


^  y  .  p  V  =  0 

0  t 


The  Naviet'Stokes  equations  for  a  Newtonian  fluid  are 

P 


ti-L 

V*  3o;J  Va.  ^)  \ 


Dv  _  3p 
^  Dt  ”  "  ^ 


h  ["  ( 


3v  2  „ 


)] 


3 

1*  i  r 

,  /  ^  +  !r  \ 

3a 

\  ^/. 

J  a.  L 

D»  _  3p  3 

^  Dt  ~  ^  ^  3tt 


i.  r  .  h.)] .  i  L  C-  *  -'ll 

3x  |_^  \3a  3yyj  BjS  ^3^  3y  ^  J 


(M) 


(1-2) 


(1-3) 


(1-4) 


With  the  assumption  of  a  perfect  gas  the  equation  of  state  ia 

p  =  p  R  Tj 


(1-5) 


where  R  is  the  gas  constant  and  Tj  the  absolute  temperature. 

Assuning  an  isothermal  film 
Tj  =  const. 

and  the  energy  equation  reduces  to 

P  a  /O  (1-6) 
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For  very  thin  filma  fully  developed  laminar  flow  exiata  in  the  bearing  gap.  Compo* 
nenta  of  velocity  normal  to  the  bearing  aurface  (w)  can  be  neglected  with  reapect  to  the 
horixontal  component  becauae  the  ratio  W/u  ia  of  the  amae  order  aa  the  angle  between  the 
bearing  and  journal  aurfacea.  Obaerving  that  variationa  of  velocity  in  the  y-direction  are 
of  order  U/S  and  U/S^,  and  are  of  larger  majptitude  than  variationa  in  other  directiona,  we 


( 


By 

Bt 


1-2,  3.  4: 

Bu 

♦  u  —  + 

=  _  11 

Ba 

Bys; 

Ba 

Bv 

v^\ 

Ie 

^  “  Bo 

ByS  / 

"  "  B/3 

_B_ 

By 


B 

+  — 


By 


Bu 

^  ’ 

By 
By  ’ 


(1-7) 


(1-8) 


_3p 

By 


(1-9) 


Ihe  left-hand  aides  of  these  equations  represent  the  acceleration  terms  which  can  generally 
be  neglected  in  conyiarisan  to  the  viscous  forces  since  their  ratio  is 


a)  Convective  Terms 

R^c  U 

^  RA)  P  U  c^  .  P  U  R 
pHr^  ■“  mR  ■  m 

c 

Ihis  ratio  is  nuch  lower  than  unity  in  most  normal  applications. 

b)  Fluctuation  Terms 

Inertia  Force  p  U'  o),  U'  c 

Viscous  Force  ~  y  p2  U  * 

M - 

c 

Using  cotimonly  net  values,  frequencies  higher  than  10®  and  full  amplitude  (u*^  -  v)  fluctua¬ 

tions  are  needed  to  make  this  ratio  of  order  unity. 


Inertia  Forces 
Viscous  Forces 


Then  the  Navier-Stokes  equations  reduce  to; 

Bp  _  /  Bu  \ 

Bo  By  \  ^  ^  /  ' 

iL  -  _L  /  Il\ 

B/3  "  ^  By  y- 


(I-IO) 


(I-ll) 


(1-12) 
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Integrating  the  e<piati<»a  of  notion,  we  get: 

Hi  -  y  y 

— *-  ■■  ■  ■■  II. 


1  Bp 


—  U; 


(M3) 


1  Bp 
"  2m  ByS 


y(y  -  hj) 


(1-14) 


From  the  equation  of  continuity,  there  results: 
B(pw)  _  Bp  B(pu)  3{pv) 


2p 

Bt 


Ba 


B/3 


Bt 


B 

Ba 


By3 


PyUj 


*‘1  HM5) 


Integration  over  the  clearnice  hj  and  use  of  the  boundary  condition  of  no  slip,  gives: 


s 


B(pw) 

Biy 


Bp 

dy  =  -  h^  + 
Bt 


i.  —  (  ^  \ 

2  Ba  \  6m  3a  / 


3/0 


/ph,3 

Bp 

\  6m 

3/0/ _ 

1  3 

2 


i.  (phu^)  -  i-L  1  (phjUj)  -  pu.  !!llI 

Ba  i  2  Ba  Ba  j 


(1-16) 


or: 


_L  1.  \  —  ( i  \  -  i 

BaXp  3a‘’/  *  BysVp  3/3/  “^j 


Bp 

2hi  ~  +  2p(wj  -  w^)  + 

y 


,  1  3h  B 

,u.  _♦  h  - 


(p[uj  +  ' 


(M7) 
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Nov,  if  n  dnotes  mgtdar  velocity: 

=  R  n. 


Wb  =  0 

3  3h, 

u.  =  Rfij - 1 

^  3t  35 

3h  j  3h  j 

Cl),  —  R  fl.  +  ■' 

J  J  3a  3T 


Then  the  right  hold  side  of  (I~17)  beconea 
6^2  h,  ^  +  2p  Rni  —  +- 


2  p 

R  n. 

3hj 

Bhi 

3t 

L 

3a 

Bt 

r 

R  n. 

3 

3hj 

Bhj 

3t 

l5  ■ 

V 

Ba 

/  ,  3  3h  \  3p 

.  K,  ^mnj  .  n,)  -  ^ 

3p  3/0  ^*>1  3 

t  h,  ma  >  Ob)^  .  2h,-  t  2p—  ^  P3- 

3  /^hi\  ^  K  i  i  !l!i 

"  ‘‘1  3^  \l5  /  3a  ■  1  ^  3a  3t  35 


(1-18) 


/^\  ^ 

\  35  /  ’3a 


(M9) 


Bit 

_L  ^ 

3t  35 

Rldj  +  (lb> 


Xj  cos  5  -  Yj  sin  5 

Wflj  +  (lb) 


(1-20) 


**1  3a  3t  35  ‘‘1  (-  »i  »in  ^  -  ti  cos  5)  _  ^ 


3  3  ®"l 


(1-21) 


C 

Then,  letting  a  =  R9,  and  neglecting  terms  of  order  — ,  we  obtain 
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m- 


=  6M(nj 


«ih«re  >  oj. 


(  2  h.)  3(p  h.)  i 

+  Cl.)  { -  - ^  +  - -  f 

**  I  (flj  fib)  3t  3^  ) 

I  can  be  replaced  by  fi 


I-A  2049-20 
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i. 


APPENDIX  II 

DERIVATION  OF  EXPRESSION  FOR  FRICTIONAL  FORCES 


Fran  Appendix  I  we  have: 

1  3p 

u  =  —  -r—  y(.y  -  h,)  + 
2/xR  3^  ^ 


(hj  -  y) 


U.  +  —  Uj 

hj  *>  hi  J 


(IM) 


Therefore, 

3u 


y  =  h, 


(II-2) 


becomes: 


BH 


cHP 


f  2R  3^ 

Then 


a  3p  ijKI  30 

- —  +  ■‘=—  +  u  — 


cH 


2v  277 

K,  =  /  SfR  cos  0  d0;  Ky  =  -  J  S^R  sin  0  d0 

O  O 

Using  the  complete  expression  for  S,  we  get: 

hi  Bp  juR^n 

K  =4  —  —  cos  0  d0  +  4  -  cos  0  d0  + 

If  t\  "iZT  •'  L 


2  30 


(II-3) 


(II- 4) 


^  -  [V  sin  6  -  \  cos  0]  cos  0  d0  -  1 1  ^2  ^  ^3  ^ 


‘1 


(II-5) 


Let  us  evaluate  the  integrals  in  (5)  one  at  the  time  using  H  -  1  -  x  sin  0  +  y  cos  0 


1  +  £  cos  (0  -  0j) 


1  1  3 

—  ^  hjp'  cos  0  d0  ~  ^  ^  ^ 

4  ph, '  cos  0  d0  +  —  4  ph  sin  0  d0  = 

2  y  1  2 

1  277  1 

—  [ph  cos  0]  1  -  —  .4  pc  [X  cos  0  -  Y  sin  0]  cos  0  d0  + 

2  o  2 


+  —  4  cp  [1  +  X  sin  0  +  Y  cos  0]  sin  0  d0  = 


Ij  =  -  4  p  sin  0  d0  -  -  4  p  [X  cos  2  0  -  Y  sin  2  0]  d0  (11-6) 
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rnm  f)  nfiR*  CO«  0  CO» 

I  =  ^2  ^  ^  - 

^  hj  e  c  coadj  1  +  c  cos  (& 


cos  6  cos 

- ; -  00  = 

.  I A  _ 


^^^2  j  cos  (9-8^)  6  sin  6  sin  6^ 

e  c  cos  8^  L'^  1  +  e  cos  (.8-9^)  ^  1  +  e  cos 

£  cos 


c  cos  8^ 


(9  -  8^) 


(8  -  8^) 

die  -8^)  + 


tan  9 


■■[‘r 


^  - - - 

1  +  e  cos  (.9  -  ) 

e  sin 


id  -  8.) 

die  -  8^)  + 


/ 


+  e  cos  id  -  9 y) 


e  cos  8  sin  9y 

id  -  8 


V 


+  ■ 

1  +  €  COS 

Bringing  the  last  term  to  the  left  hand  side,  we  have; 

RnP,  cos0  "1  , 

- I  — 1 - rs - 2~v  \  ^1^ 

c  I  1  +*  €  COS  (6  ^  6^)  \  * 


/uR 


€  C  COS  6y  I 


€  COS  id  -  9  y) 

- - - -  & 


id  -  dy) 


-  tan 

Tliere  fore: 

I2  =  — - 

‘  ec 


^  - 

1  +  f  c  o  s 

ie  -  dy) 

die  -  By) 


"]■ 


€  sin 


e  ^  - 

1  +  e  cos 


ie  -  8y) 


cos  9  y  d0  -  ^ 


d9 


+  e  cos  id  -  9y) 


tan  9y  4  d  In  [1 
27T  "I 

ry 


+  e  cos  id  -  6 


VI 


-  COS  a, 

€C  * 


I  277 


tan 


I 

2  e  c 


L  rr~jj 

r  _  —  -4 . -1y 

e  c  |_  n  “ 


277  > 

9  [In  [1  +  e  cos  id  -  8y)]  ]  | 
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,  _  •>  fA  ^  ^ 

I  •  ~  IV^Y  — 66  ~ 

1  +  «  coa  (6  -  dj) 

=  u —  I  4  -  d5  = 

2  ■  1  +  £  coa  (5  -  5j) 

=  M  X  »i"  2  Oi)  ^ 

2  '  1  +  €  cos  ^ 

where 

^  =  5  -  5i 


f/RV  1  /«/vo  9 

5, 

sin  i  cos  ^ 

HA  •f 

fini  1  cos  z 

1  +  c  cos  f 

05 

sin  2  6. 

2  ^ 

cos^  f 

sin^  ^ 

d.l 

1  + 

£  cos  ^ 

d# 

27T  -  6^  ^  ^ 

/  ^ - L  =  0 

-5j  1  +  c  cos  g 

because  if 
sin  ^  cos  i 


«ASB  ^  ^^V/O  V, 

“Tcos^  = 

fj  (f)  =  -  (-f)  (odd) 

Then: 

I3  =  M  -r  ^  sin  2  5.  r  .4 


_d£ _ 

e  cos  f 


-  2  4  - 

^  1  + 


lin  2  6t 


e  cos  ^ 


{  i 


R  .  n  ,  sin' 

^2  ^  ^  1  +  e  cos  f  1  +  e 

=  ^  -i"  2  «1  [♦  i  .  i  ,♦  If 

-  i  *  o..f  «  .  _J£_|J= 

R  .  f  /  2  -  eV  277  477  “1 

= "  i  -  ttJ 


sin^  ^ 


e  cos  ^ 


/X  —  t  sin  2  5 

2 


1  -  e 
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-HTirV 


M  R"  2XY 


£*  ■!  1  - 


[^I 


(ii-e) 


I.  =  -  — 


COS 


e  cos  (6  -  5j) 
the  last  term  to  the  left  hand  side,  we  get: 


.  ,  cos  2  8  +  1  -  cos*  8 
68  =  -  /xRJ^  -r- - rr: - r-r -  68 


1  +  e  cos  (6  -  8j) 


21^  =  -  /xRX  ^  — 


cos  2  0+1 


€  cos  (8  -  0 j ) 


d0  = 


=  -  ;xR)(  L  ^  ^  _^L_1 

1^^  1  +  e  cos  ^  1  +  £  cos  if  J 


xRJC 


[• 


2  0,  4  - 

1  +  1 


cos^  ^  -  sin^  ^ 


€  cos  ^ 


6f  + 


,  sin  i  cos  (f  , 

-  2  sin  (20.)  4i - ^ ^  d^  -  mRJ( 

*  1  +  €  cos  ^ 


-J 1  -  e* 


=  -  /xRX 


l^  = 


P-;"?'*  (i-jrr;^)^  .  ...  ..> 

L  €*  J 1  -  €*  J 1  -  e*  J 

'i  _  n~r7^v  Y*  -  x* 


JTT7? 


T 


+  1 


Gmcluding; 


(11-9) 


K  =  —  {  R  P.  4  P  sin  0  d0  +  R  P  4  P(Y  sin  2  0  -  X  cos  2  0)  d0  + 
*  2R  J  •  • 


AR  P. 


fiTT* 


/  c\  I  X  /  1  -  H  -  V  /  2  XY  y  -  Y*  i  +  X*  i\  I  1 

U)|-«'  (— ( - 77 - )u| 


(II-IO) 


Now,  for  ky 


-  4  —  ““  sin  0  d0  +  4  — -  sin  0  d0  +  -  uRC  4 

2  30  ^1 

sin  0  d0  =  Ij'  +  I2'  +  I3'  +  I4' 


[ 


t  sin  0  -  X  cos  0 


] 


(IMl) 
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(IM2) 


h'  = 


ttA  I'P*  ( 1  -  J  1  -  c*) 


”  ^  rrrp- 

mRY  j^(Y2  -  x2)  ^ 


1  -vTTTT? 


)  -  ‘0 


^  <1  -rrr?-) 


(11-13) 


(II-U) 


(IMS) 


J  1  - 


Therefore: 


K„  =  <R  P.  M  coo  a  d0  +  R  P.  >  P  [Y  cos  2  0  +  X  sin2  i9]  de  + 

y  2R  I  ■  • 


ir 

+  — 

3 


J  1  - 


(11-16) 


Considering  that <  <  1.  ~  0(1).  ^  ~  0(1)  the  expressions  for  and  k,  become: 


Kx  = 


"  (^)  1 


277 


^  P  sin  0  60  *  P  [Y  sin  2  0-  X  cos  2  d^  + 


2t7  ^  (1  -  h  -  e2) 


°  (i)  1 


e2  ^1  -  e2  I 

^  P  cos  0  60  +  P  [Y  cos  2  d  +  X  sin  2  0]  d0  + 


(lun) 


2it  ( 1  -  '1 1  -  e^) 

+  -  A  - -  X 

^  £2  frr;^ 


(IM8) 
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APPENDIX  III 

SOLUTION  OF  TRIDIAGONAL  SET  OF  N  SYSTEMS  OF  TWO  EQUATIONS 


Consider  the  system: 

“i  *i.l  +  ‘‘i  *i  *i+l  *  ^  *i 

h  Vi-l  *  *  *'i  Vi  +  l  ^  ji 

i  =  1  -  N 

Since  the  terms  x^,  are  kno«n,  in  genersl  «e  hsve  relations  of  the  type: 


(III-l) 


>1-1  '  *1  -1  ♦  ,,  t  C| 


'  I'l-l  >  “i  <i  *  Bi  ’l  * 

Then,  by  substitution: 

Xi(bi  +  A. a.)  +  ^i*i  +  l  "  *i  ■  ^i“i 

yi(«i  +  E.fi)  +  hjy.^l  +  x.{ii  +  f^D.)  =  ji  -  F.fi 

These  equations  can  be  expressed  as 
X.  W.  +  yj  Z.  =  Ri  -  cj  xj  +  i 

*i  ^  yi  ^  -  ‘‘i  yi  +  1 

Or: 


Ri  -  cj  Xj  +  i  Z. 

W.  R.  -  Ci  Xi^i 

Si  -  ^  yifi  Yi 

Xi  S,  -  h.  y.,1 

w.  z. 

’  Yi 

Wi  Zi 

Xi  Yi 

Xi  h 

[herefore: 


YjRj  -  Y,c,x.^^  -  Z,S,  ^  Z.h.y,^, 
WiYi  -  X^Z. 


(111-2) 


(III-3) 


(III-4) 


-  W.Sj  -  W.h.y^^i  -  RjX^  ^ 

Vi  - 


(111-5) 
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Comparing  with  III-2,  we  have 


i+l 

=  -Y.c./Lii 

Di+l 

*■  X.c./l. 

=  Z.hj/Li  J 

Ei+i 

=  -W.hi/Li 

i+l 

=  (RiY.  -  Z. 

S,)/L,; 

Fi+1  = 

«4iere 


(III-6) 


X. 


S,  = 

1 


+ 

‘‘i 

+ 

•iBi 

e. 

1 

- 

«iCi 

8i 

+ 

f.E. 

ii 

+ 

Ji 

- 

WjY 

i  ~ 

x.z 

(111-7) 

After  evaluation  of  A.,  Bj,  C.,  Dj,  Ej,  Fj,  for  i  =  1  -*  N  +  1,  use  of  (111-2)  for  i  =  (N  +  1) 
”*  1  will  give  the  answer  to  the  problem. 
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Ihe  prograaned  equations  were: 

=  X 

4U.R,  2K, 

-  B  ^  P  sin  ^  d5  +  -  cos  fiT  + 

ss  .  /^s 


=  0 


B  P  cos  0  dd  +  BL 


Me 


sin  nr 


« -  *.l, 


(■ 


.X  + 


[Y 


+ 


Kj  and  are  the  artificial  damping  coefficients  in  the  X  and  Y  directions  respectively, 
a^  and  major  axes  of  assumed  elliptical  idiirl. 
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C  INF.  JOUR.  8RG.  V.  CASTELLI  FRANKLIN  INST. 

10  DIMENSION  WllOOIf  WitlOOIt  HllOOIt  0TI20)t  STHETAIlOO) 

DIMENSION  CTHETAMOO),  OHDTHEIIOOIt  W2(100) 

21  DIMENSION  OXDT1M20).DYOT1(120).X1(120)»Y1(120)fNSTn20) 

DIMENSION  X2(60001.Y2(6000).LP(S1»51 1 

FREQUENCY  90(101.3401  301.4201  1.101.4601  1.  10 1 . S 1 0( 30. 1 . 0 1 . S50( 30) » 

16  30(30)  .660(0.1 .100.670(0. 1.  n.730(  100. 1 . 0 ) .  740( 0. 0. 1 )  .763(30). 
1810(20. 1.0) .830 ( 1. 100). 900( 1. 1,0), 920 ( 1.1.0) 

1  CALL  LGCHAR(S,4HUX03) 

30  RE AC  40,N,NOT,L1.L2.L3,L4,LNSTEP,LKOUNT,NPO.NSPA,TRUNC.DXOTO.OYOTO 
I.NCASE 

40  FORMAT (213, 12, 3 1 1 , 14, 1 3, IS, 1 2, 3E 14.8. IS) 

50  F.N  *  N 

60  OTHETA  »  3.  14159265*2. 0/FN 
70  THETA  ■  OTHETA 
80  Nl  >  N  *  I 
90  DC  120  I  >  2,  N1 
100  STHETA(l)  *  SINF(THETA) 
no  CTHETAd)  *  COSF(THETA) 

120  THETA  =  THETA  ♦  DTHETA 

121  DTHE2  *  DTHET/.*OTHETA 

122  D2THE  *  2.0*0THETA 
1.30  SENSE  LIGHT  0 

140  1F(L1)  (60,  (60.  130 

150  SENSE  LIGHT  1 

160  IF(L2)  180,  180,  ) 70 

170  SENSE  LIGHT  2 

180  IF(L()  200,  200,  190 

190  SENSE  LIGHT  3 

2C0  IF(L4)  220,  220,  2)0 

2  10  SENSE  LIGHT  4 

220  READ  2.30,  (DT(K),  K  =  1,  NDT ) 

230  FORMAT  (5E14.8) 

240  REAU  230,  PLAMOA,  B,  PL,  XO,  YO,  UNBAL,  DAMP1,  0AMP2,  ABRTIO.BNUMB  . 
250  READ  230,  (W(I),  I  =  2,  NI) 

C  PRELIMINARY  CALCULATIONS  OVER.  NOW  SET  VARIABLES 

260  T  =  n.o 

261  L=C 

262  LL=0 
270  X  =  XO 
280  Y  =  YO 

290  K  =  1 

291  ASSIGN  661  TO  Ml 

292  NSTEP  *  0 

293  ASSIGN  525  TO  M2 

294  KOUNT  =  0 

295  8N  =  0.0 
300  DXDT  *  DXDTO 

310  PYin  =  OYOTO 

311  BPL  «  B  *  PL 

C  END  OF  SETTING.  ENTER  MAJOR  LOOP 
320  SUMl  aO.O 
330  S0M2  *  0.0 
340  DC  .390  I  «  2,  NI 

350  H(I)  =>  1.0  ♦  X^STHETAII)  ♦  Y*CTHETA(I) 

360  OHDTHEd)  =•  X*CTHETA(I)  -Y*STHETA(I) 

370  p  »  w( I )/H( n 

3fe0  SLM1  *  SUMl  ♦  P*SThETA(I) 

390  S0M2  =  SUM2  ♦  P*CTHETA{ 15  A 
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400  G  >  SUH1«0rHETA*B 

410  0  >  SUM2«0THETA*A  ♦BPL 

420  IFISENSE  LIGHT  1)  430*  460 

430  SENSE  LIGHT  I 

440  G  «  G  «  UNBAL*C0SFI2.0*T) 

450  Q  >  Q  >  UNRAL*SINF(2.0»T) 

460  IF (SENSE  LIGHT  21  470.  500 
470  SENSE  light  2 

480  G  >  G  ♦  0AMP1«(ABATI0»(Y  -  YO)  -  OXOT) 

490  Q  «  0  -  0AHP2*((X  -  XO)/ABftTIO  *  OYDT) 

SCO  T  -  T  ♦  DT{K) 

510  IF(r  -  97.3M93722)  521,  520,  520 

520  1  »  T  -  97.3b93722 

521  DXOT  »  OXOT  ♦  G«OT(K) 

522  OYOT  »  OYDT  ♦  0*DT(K) 

523  X  »  X  ♦  DXOT*OT(K) 

524  Y  »  Y  ♦  OYOT*OT(K) 

5241  L  *  L  ♦  1 
LL  «  LL+1 

5242  OXOTl(L)  »  OXOT 

5243  OVOTKL)  «  OYOT 

5244  XKL)  a  X 
X2(LL)»X 

5245  YKL)  »  Y 
Y21LL)*Y 

NST(L)  «  NSTEP  ♦  1 

5246  IF(L-100)525. 5247, 5247 

5247  WRITE  OUTPUT  TAPE  15, 524B,N<CASE ,  b,  (NSTI L  ) ,  X  |  (  L  1 ,  Y 1  ( L )  ,  OXOT  1  (L 1  ,DYD 
ITKLl.NSTIL^SOl.XKL^SOl.YKL  +  SOl.DXOTKL+SOl.PYDTKL^SOl.L-l.SO) 

524b  FORMAT! 12HICASE  NUMHER  15, 4H,  B«  ElS.d/6H0  STEP  SX  IHX  lOX  IHY  lOX 
15H0X/0T  10X  5H0Y/0T  IIX  5HSTEP  5X  IHX  lOX  IMY  lOX  5H0X/0T  lOX 
I5hOY/OT//(I6,2F11.7,2E15.7,5H  **  15,2F 1 1 . 7,2E 15.7) ) 

5249  L  =  0 

GO  TO  525 

525  SUMS  »  0.0 

526  SUM4  «  0.0 
530  M(1)  *  WIND 
540  H(N  ♦  2)  «  H(2) 

550  00  620  I  »  2,  N1 

560  GAMMA  >  IMd^n  *  MII-1)  -  2.0«MI  I )  ) /0THE2 
570  CSI  =  IH(I*n  -  W(  I-1))/02THE 

580  FOT  «(-(H(I)*(W(n  +  CSI  *  OHUTHEI 1) )  -  H( I )* (CSl*CSI  ♦  WII)»( 
1GAMMA  ♦  W(n)))/PLAMOA  -  CSI)*DT(K) 

60C  SUM3  «  SUM3  ♦  ABSF(FOT) 

610  WKI)  -  Wd)  ♦  FOT 
620  StM4  »  SUM4  ♦  WKI) 

629  NSTEP  »  NSTEP  ♦  1 

630  00  640  I  -  2,  N1 
640  W(I)  «  WKI) 

650  CCNINO  «  SUM3/SUM4 

660  GO  TO  Ml,  (661,  665,  SCO) 

661  ASSIGN  665  TO  MI 

662  GC  TO  730 

665  WRITE  OUTPUT  TAPE  16, 666, NCASE, NSTEP, CONINO 

666  F0RMAT(9H  CASE  NO.  15, 9H, STEP  NO.  15, 5H,  M  »  E15.8) 

670  IFICONINO  -  PRECM)  730,  680,  680 

680  IFISENSE  LIGHT  3)  710,  690 

690  SENSE  LIGHT  3  A 

700  GO  TO  730 
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710  IF(K  -  NOT)  720.1320,I320 
720  K  ■  K  ♦  1 

710  IKNSTEP  -  LNSTEP)  7%0t750*  750 
7>i0  1F(C0NIND  -  T%UNC)  750.  750*  770 
750  ASSIGN  800  TO  Hi 
740  ASSIGN  890  TO  H2 

761  WHITE  OUTPUT  TAPE  99.  762.  NSTEP*  NSTEP 

762  F0KMAr(26H0  START  OF  ORBIT*  STEP  NU.  IS/I0H2  IN  ORBIT  IS) 

763  00  764  I  -  2.  N1 

764  W2(I)  «  Nil) 

765  GO  TO  820 

770  PRECM  «  CONINO 

800  KCUNT  «  KOUNT  *  1 

810  IF (KOUNT  -  LKOUNT)  850.  820.  820 

820  KOUNT  >  0 

830  WRITE  OUTPUT  TAPE  I 6.e40.NCASE.NSTEP . B* T.X. Y.OTIK) . PLAHOA, I W( I ) » !• 

12. N1) 

840  F0:tMAT(6H0CASE  14. 6H,  STEP  I5.4H.  B>  E14.7.7H.  TIME-E14.8(4H»  X« 
lr'll.8.4H.  Y«F11.8.  5H«  DT>  E14.8/9H  LAMBOA>E  14.8*  10X  22H  W  >  PH 
1  DISTRIBUTION/! 1H  I0FI1.7)) 

850  GO  TO  H2. (525.890) 

890  EPSIL2  «  X»X  ♦  Y»Y 

900  IF(6PSIL2  -  0.98)  910.  9)0.  11)0 

910  BN  a  BN  4'  1.0 

920  If (8N-BNUHB)320.930.930 

930  BN  a  0.0 

931  IF(SENSE  LIGHT  3)  9310*9310 

9310  IF(SENSE  SWITCH  2)9398.9311 

9311  BIGX»X2(1) 

9312  DC  9315  LLla2,LL 

9313  IFIBIGX  -X2(LL1))9314, 9315. 9315 

9314  8IGX  a  X2(LL1) 

9315  CONTINUE 

9316  SHALLX»X2(1) 

9317  DO  9320  111-2. LL 

9318  IF(X2(LL1)-SHALLX)  9319,9320.9320 

9319  SHALLX  «  X2(LL1) 

932O*C0NriNUE 

9321  BIGY-Y2(1) 

9322  DO  9325  LLla  2,LL 

9323  IF1RIGY-Y2(LL1) 19324,9325.9325 

9324  OIGY  a  Y2(LL1) 

9325  CONTINUE 

9326  SHALLY  aY2(1) 

9327  00  9530  LL1a2,LL 

9328  IF(Y2(LL1)  -  SHALLY)9329, 9330, 9330 

9329  SHALLY  •  Y2(LLl) 

9330  CONTINUE 

9331  IF(SHALLX)  9332,9333,9333 

9332  SHALLX  -SHALLX  -  0.1 

9333  LOWX  a  SHALLX  •  10.0 

IF (SMALLY)9334, 9335, 9335 

9334  SHALLY«SHALLY-0.1 

9335  LOWY  a  SMALLY»10.0 

9336  SHALLX  •  LOWX 

9337  SHALLY-LOWY 

9338  SHALLX  -SHALLX/IO.O 

9339  SHALLY  ■  SHALLY/10.0  ^  ZQ 

9340  LANGEX-  I RIGX-SHALLX)/0.2  4-1.0 
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L*NGEV»  (BIGY-SHALLY)/0.20«1.0 
IF JLANGEX-LANGEY 19343,9344, 9444 

9343  RANGE  -  LANGEV 
LANGE-LANGEY 
GC  TO  9345 

9344  RANGE>LANGEX 
LANGE«LANGEX 

9445  RANGEY«0.20*RANGE 
RANGEX«0.20«RANGE 
eiGX>SMALLX«RANGEX 
BIGY=SMALLY*R4nGEY 

9347  NPLO=ILL*NPO-n/NPO 


9348  LLIM«''iPLO»,'IPO 

9349  DC  9451  I»LL,LL1M 

9350  X2(l)  =  X2{LL) 

9351  Y2(n  =•  Y?(LL) 

DO  9362  NPL*1,NPLO  dH  /) 

CALL  FNPLOr  (29H(24H  INF.  BEAR.  SHAFT  ORBIT)  ,  8H(--3H — X) 
lV») ,SMALLX,BIGX,SMALLY,eiG¥,20, 1 ,20, 1 ,6Ht F 5.2 ) , 6H( F5.2 ) ) 
L2L-NPL*NPO 


3M. 

,8H(  3H 


XXXXX-SPACEF(NSPA) 

CALL  CURVE(X2(L2L),Y2(L2L),NPO,6H  ) 

WRITE  OUTPUT  TAPE  99,9361 ,NCASE,B,NPL 
9461  F0RMAT(9H  CASE  NO.  14, 4H,  B»  E 15.8, 9H, PLOT  HO.  14) 

9362  CONTINUE 

9398  LL»0 

9399  BN»0.0 

940  read  950  ,  B,  BNUMH,  INDEX,  LI,  L2,  L4 
950  FURHAT(2E14.8,  15,  312) 

960  IFISENSE  LIGHT  1)961,961 

961  IFIS6NSE  LIGHT  2)962,962 

962  IFtSENSE  LIGHT  4)970,970 
970  IF(Ll)  990,  990,  980 
980  SENSE  LIGHT  1 

990  IF(L2)  1010,  1010,  1000 
1000  SENSE  LIGHT  2 
1010  IF(L4)  1021,1021,1020 

1020  SENSE  LIGHT  4 

1021  WRITE  OUTPUT  TAPE  15,  1023,  INDEX,  L),L2,  L4 

1022  write  output  TAPE  99,  1023,  INDEX,  Ll,L2,  L4 

1023  F0RMAT(9H  INDEX  »  14,  312) 

1024  BPL  =  B«PL 

1030  GO  T0(  1031,  1031 , 1091, 1031  ), INDEX 

1031  IF1L)1033, 1043, 1032 

1032  WRITE  OUTPUT  TAPE  15, 5248,.NCASE,  8,  I NSTI  IN)  ,  X 1 11 N) ,  Y  1 1  IN)  ,  OXDT  1 1  IN) 
1,DY0T1 ( IN),  IN  *  1,  L) 


L»0 

1043  GO  TOI 1040,40,1091, 1260), INDEX 
1040  DO  1050  I  =  2,  Nl 
)050  Win  *  W2(  I  ) 

1060  X  3  XO 
1070  Y  =  YO 
1080  DXDT  -  OXOTO 

1090  OYOT  »  DYCTO 

1091  IFISENSE  LIGHT  4)1092,1100 

1092  READ  40,N,NOT,H,L2,L3,L4,LNST£P,LKOUNT,NPO,NSPA,TRUNC,OXOTO, 
lOYOTOiNCASE 

1100  IFISENSE  LIGHT  3)1340,320 

C  failure  PROCEDURE. SIGNAL  ON  THE  PRINTER  AND  READ  NEW  INPUT  DATA. 
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ino  WRITE  OUTPUT  TAPE  l!}.  1 1 30*NCASE •  B*  NSTEP 
1120  WRITE  OUTPUT  TAPE  99,n50t  NCASE.  B»  NSTEP 

1130  FURHAT(26H  BEARING  FAILURE.  CASE  NO.ISfWH,  R-  E 14. T|9H.SrEP  NO. 
1IS/IH2I 

1140  DO  1160  I  •  2.  N1 
IISO  Win  >  W2(I) 

1160  X  >  XO 
H70  Y  -  YO 

n«iO  Dxur  »  oxDTc 

1190  DYDT  ■  OYOTO 
lifCP  GC  To  '’^1 

C  TE*?HIM*TI0N  PROr^DURE  TOR  END  OF  JOB.  REWIND  AND  PAUSE. 

126o  »JRI(E  OUTPUT  TAPE  99.  1270 

1270  FURNAT(29H  JOB  FINISHED,  MACHINE  PAUSE.) 

12SO  ElvD  file  lb 
l?:i  END  FILE  16 
l:-'”'’  CALL  LGCHAfUS.'tHUXCi) 

FNf  FILE  5 

1290  REWIND  IS 

1291  REWIND  16 
1^92  REWIND  b 
l^rn  PAII^E 

IblC  CALL  E.,DJOC 

C  PROCEDURE  FOR  CASE  WHEN  LIST  OF  DELTA  T  IS  TERMINATED. RESTART. 
1320  WRITE  OUTPUT  TAPE  99,  1330 

1330  FORMAT (45H-EXHAUSTEC  LIST  OF  DELTA  T.LOOK  FOR  NEW  CASE./1H2) 
Ij40  sense  LIGHT  3 
13S0  GO  TO  9U0 
END 
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X/OFiwiTt.  )oo 

iFjfJi  7€  Ojifcir  l^yrofiSfi-M 


Ol  oSS/lW 


A/"  *  no.  oi^  'poto^’s.  «iroO'<i4 
W^T  =  /9o.  ztApid 
LMSfep=  ^n 

LIcou^st  »  7>r>W^  'pressure  cij:fer  i-icoot<fT  s^i^s 

slpo  JT  nu»Ji.«r  of  X-V  p-Ai’rs  p'tcTa.rc 


0  ' 

AJS?fl=  -p^Df  €.U«r^  AJiTft  p>oi'»J't- 

yi^OKic  ■=■  IK  c>t-t>i'f'  uikeK.  IS  ^tss'f^K 

pKDTo  r 

JPVDTo  =  ^ 

Flam6»=  a 


B  -  S 

XO  =  X„ 

Vo  -  t:  . 

yUBftu  =  ini' 

M  C 

"bfiM?!  =  '2- k;/(^X2.M} 

T>AH?2  ^  2 

/16^r(0  *  f  m't’o  o^  ciXfS  eiSSeJ^MtJ  edJUfu. 

3n4/W?>  ■=  f^.oj  Tstcf>S  Le^ort  neKT 
IV  -  Ty  /)Z3 
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Iteluia 


««  S/U^ 


ChUf  of  UnX  ItMank 
Daia  rtaiiit  af  %)ia  iaagr 
Vaahli«|«ii  tS.  0.0. 

Attni  Co4a  4M  5 

Attat  Coia  Mt  1 

Mtai  CoAa  kil  \ 

Attnt  Ooda  09  1 

Aiini  CoAa  AM  1 

rn—aiirtlm  Offiaar 
Offlaa  of  lawa  Kaaatrah 
iyaiMli  (tffioa 
A9S  SwMT  Siraai 

iootoA  10,  Ihaaaohnaatta  1 


CoMaadlAg  Offloor 
Offteo  of  laml  laaaaroh 
iNBeh  Offieo 
06  laat  KaaAelph  Straot 

Chleage  1,  niinola  1 


CoMaadiai  Offloor 
Offieo  of  lawl  laaoaroh 
Braiwh  Offieo 
307  Moot  34U  Strooi 

Mom  York  U,  M.T.  1 


CaaaanAlfig  Offloor 
Offieo  of  Maml  loooareb 
■roneh  Offioo 
1000  Ooory  Stroot 

Son  Froneiooo  9,  California  1 


CflMonding  Offloor 
Offieo  of  Haoal  laaoaroh 
Brofieh  Offieo 
1030  Boot  Oroon  Stroot 

Pooodono  1,  California  1 


Chlof .  Buraaa  af  laial  Vaapina 

DopartaiwA  of  tiia  lla«r 
Iteahlngton  Mi,  0*0, 

Atint  OMaa  MMt  (J*  V*  Oaaa) 


1 

I 


Diraotor 

Natal  Roaoaroh  Ubwalary 
Matitinfton  IS,  0*0* 

Attni  Coda  1000 
SI90 

Spoeial  BroSaota  Offiaa 
Dopartnent  of  tha  Baop 
Voahlngton  IS,  D*C* 

Attni  Coda  SK3*4  (0.  Gold) 


6 

1 


1 


Head,  Boorlnga  and  Stala  Iraneb 
U.S.  Matal  Inglnooring  fc^irtnint 
Station 

Annapolio,  Nairlaad 
Attni  Coda  ISl  (Vatt  T*  Snith) 
Ifr.  Bcbart  BoU 


1 

1 


Notorial  Laboratory  Lttaarjr 
Building  191,  Coda  9ia 
Bom  York  Itoml  Bhippard 

Brodtlyn  1,  Vom  York  1 


Library 

Toehnieol  Xoporte  Saetien 
U.8.  Nfttol  Boatgroduoto  Sehool 

NontM^,  Collfemla  1 

Ceonandiag  Offloor 

U.S.  Naval  ArionieB  Paeility 

ZndlaMpolla  U,  Indiana 

Attni  J.O.  Voir  1 


Ceanandlng  Offlear 
Offieo  of  Naval  loooaroh 
Broneh  Offieo 
Ben  39.  Navy  #100 
Bloat  Beet  Office 

Nom  York.  Mom  York  10 

Chief,  Bureau  of  Shipo 
Oepartaent  of  the  Mavy 
Maahingten  35.  O.C. 

Attni  Cede  6AAB  (Jaaao  C.  Bold.  Jr.)  I 

Attni  Code  343  (I.  A.  Bukain)  1 

Attni  Code  634  (B«  A.  Vapaa)  1 


Direeter 

U.S.  Naval  Bailor  A  Turbine  lab. 

Naval  Bare 

Philadelphia  12.  Bwiaylvaaia  1 

Offieo  of  Chief  of  Ordnaneo 
Beaeareh  and  Dovelepant  Dlv. 

Oopt.  ct  the  Any 
Vaohlngten  25.  D.C. 

Attni  Meraan  L  Klein  1 
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PueU  k  Lubrieante  Section 
Besearch  Branch 

Reeeareh  A  OorelepMnt  Oivioion 
Office,  Chief  of  Ordnaneo 
4th  Floor,  Bentagon  Aimmk  #2 
Waohington  25,  O.C. 

Attni  Ik*.  Ronald  B.  Btreoto 

Chief  of  Roeoareh  and  Develepaeiit 
Office,  Chief  of  SUff 
Oepartaent  of  the  Amy 
Pentagon  Building 
Vaohlngten  25.  D.C. 

Cceanvilng  Ooneial 
U.S.  Antr  Bnglnoor  RAD  Laboratorioe 
Bert  Bel voir,  Virginia 
Attni  V.H.  Crln,  Nvieloar  Bawpr 
Field  Offfioe 

Director  U.S.  Amy  Bngineer  Botoareh 
and  Developaent  Laboratory 
Fort  Bolvoir.  Virginia 
Attn:  Teehnieal  Deemonto  Cantor 

Comander 

Amy  Rocket  A  Guided  NiooUo  Agamy 
Redaiene  Areonal.  Alabom 
Attni  Taehnieal  Library 

BSD  -  AROD 

Bm  CM.  Duke  SUtion 

Ourhaa.  North  Carolina 

Chiof  of  SUff.  U.B.  Air  Foroo 
Tho  Bmtagon 
Waahlngton  IS.  D.C. 

Attni  ABBBUAB^ 

Comandor 

Air  roNO  Offieo  of  Seiontifto  leooaroh 
Waohington  15,  D.C. 

AtUi  SNMI 


Amod  Sorvleoo  Toehnlool  Xnfcraatlm  Agemy 
Arlington  Hall  SUtien 

Arlington  12.  Virginia  10 

Offiee  of  Aooiotant  Direeter 
(Aiegr  Roaetoro) 

Dlvlolon  of  Reactor  Deaelepamt 
U.S.  Atenie  biorgy  jCemiooien 
Vaehlngten  25.  D.C. 

Attni  Nr.  Claromo  I.  NiUor.  Jr.  2 

Nr.  N.  Oroeonan,  Chief 
BnglnoerlAg  DoTolofmnt  Braneh 
Reactor  Dovolopmnt  Divioion 
U.S.  Atenie  btwgy  Cemiooion 

Vaehlngten  25.  0.0.  1 

Hoadquartera  Library 

U.S.  Atoale  Inorgy  Ceaniaoion 

Vaohlngten  25.  D.C.  1 

Chiof.  Divlaien  of  Inginooring 
Naritiao  Adainlotration 
QAO  Building 

Vaahlngten  25.  D.C.  1 

Ciyegenio  Bnglnoorlng  laboratory 
National  Bureau  of  Standardo 
Boulder,  Colorado 

Attni  Library  1 

Rr.  Harold  Noooing 

Natienal  Aorenautioo  and  Bpaoo  AAHaiotratioa 
1512  H.  Stroot,  N.V. 

Voohington  25.  D.C,  t 


Nr.  Idaund  B.  Biooon 

Chiof.  Lttbrieation  A  Wear  Bimneh 

LiMio  Raaoaroh  Cantar 

Batienal  Aorenautioo  and  Bpaoo  Adilaidtiatioa 
21000  Broekpark  load 
Clovaland.  Ohio 


1 


Ccaaandor 

Aoronautieal  Byrtoa*  Divlaior. 
of  the  Air  Boreo  lyiftoao  CeMand 
Wri#it-Bbtteroon  A#  Baoa,  Ohio 
-  ‘  A8II0NL-2.  J.L.  Ilwrid 

- ,  B.D.  Naalon 

•1.  O.L.  NoBaddm,  Jr. 

>-i,  i.v.  1 

,  0*A*  Oaana 


lb>.  RudoIMPh  Bayar 

Oaorgo  C.  Nudhall  Spaoo  Blight  Coatar 
V  National  Aoronautioa  and  Bpaoa 

Adainiotfation 

Ouidanea  and  Control  DlvSdion 
QyrofcSUbilieer  Oranoh 

Kunumio,  Alibam  U 

Nr.  H.V.  lavage 
Oik  Ridgo  Natioml  Ubemioiy 
Boat  Wioo  l«K  Y 
Oik  Ridgo,  Toanoai 


3 


4 


1 


r.P.  tall 

QlkMta. 

tatUi 


Mfti.  ta 


ilpUal  flMlM  Utmm/m 

ISmUtfim  MMMttr 

tita  lirtit*  Ifciilliil 


Bipfta*  if  OtalMl 
taTtfktannlfer 
ta  T«Mt  9*  ta  M 
Mint  Am  J.  liilMrt  iMt, 

«r  iMtar  tatartai 


PrrfMtr  A.  (taM 
fiM  VwMliflMl 

taMiM.  lUiMdi 


X 


X 


X 


1 


MUM  taMMAM 
orriM  tfAmikMMh 
•/•  oaiMAlAy  9t  tantaali 
34)l«UiMAMrM 

mtalfM  4.  tai^  4 


FrrfMtr  F.  &• 
ToHM  SMAX 


ttaiMVity  Nm. 

FhiXtalAAU.  PMi«  1 


Jai  FramXiiw  trtMtnw 
OalifamU  MiUoU  *t  ImMUv 
AMO  QU  ITM  M. 

FMiM.  OalifMU 

Aitat  Hr.  tarl  Aim  X 

AtiAi  UbiMy  X 

IMUM 

Chlaat*  nUMt 

Attai  FrrfMit  t«l.  Ta«  X 


Mmot  11.0.  SM,  ta« 

^  Hwliaidiit  latfniwim 
Cmv^U  batiWU  9i  TwhartiT 
FitUAor^  U,  taMflMU 


taUMCiac  FNjMU  UAwWiT 
NaMMlMMita  iMiltaU  «t  tMfeMUff 


■atUUA  WmmitX  tailM 


Attoi  Or.  ImamU  O^ta 

taklU  fiwUlMlr 


Utatarr  fir  Ummk  mX  BMltai* 
PtiilMlgMi  “ - “ — 


Attfti  Frafitiir 


taufli 

mtOId. 


foUar 


library  ZnatitaU  af  AaraapMA  laiMM 
2  laat  lAAh  SkraaA 


Mr.  O.B.  SyaM) 

Ar.  Haabar.  fattalaaX  AUff 
nf  Falaral  Ubarabarlaa 
Olrtaiai  af  taMbiaMl  Tatafbaaa 
aaA  taM^  OaryaaaAUi 
151S1  IXairaa  AiraaA 

,  OaUfamU 


Aarojat  QanaraA  Mwlawlaa 
y.O.  lev  il 
Aan  Amr*  OallfMla 
Atiai  Barbara  N.  Frabart 


Mraapaaa  CaryaaaaiaR 
F.O.  B«r  9S0B$ 
laa  AS.  Oalifaarla 

Attoi  Aarearaaa  Ubaary 
Aitat  TaahAiaal  Baparta  Ormp 


AlAaaaarah  MHMfbatarliM  OMaay 
M  Marbar  AlrBaat 
AOA  Aattib  )6tb  BtrM 
fboMa.  Ariaaaa 
atiai  LibaarUa 


Mr.  VlUlaa  D.  AtMal 
Aaaaarah  Ubaratarlaa  Library 
mia  rtialaaaa  MMlbatarlat  Oa. 
Milwaabaa  1.  ViaaaMia 


J.W.  Faabar  Divlaian 
Aaariaan  CMiaaX  Oa. 

A709  Bata  BXU. 

FltUburfh  9.  FiBRAfl«ala 


•  S  • 


1  Amrleaa  SooiaV  af  babrlaatUa 

Sa(iaa«ra 

$  Barth  Vabaab  Aaaaaa 
Chleafo  2,  UXiaaU 

2  ChalMB 

Raaaarch  CaaaitUa  an  Ubriaatlaa 

Tha  Mriaaa  AaiUtf  af  IhabiBlail  tailaaara 

Unitad  BagiMarUg  Oaatar 

3l$  Saat  ATtb  Atraat 

Naw  Tork  IT*  Baa  Taib 

I  Nr.  Jaaaa  A,  Karr,  MaUaai 
l^nMRlng  OiTltlaa 
AVOO 

S*ratfard(  CatataeilaaB 

1  AaaaarQh  Fraeiaieo  KaohaaiaM 

DiviBion  of  tha  Bardaa  Oorparatlaa 
A  Old  Naatowa  fbad 
Dantmiy,  Cacnaatlaat 
AttBi  B.L.  NUi|  Tiaa  FraaUaoA 
Oanaral  Naaaiar 

1  Bamar  Enginaarlag 
Xnduatrlal  Fnit 
Fart  WaahliMtaa,  PaanayVaaaXa 

1  Utiea  i)iviaian 

Tha  BandU  Carparatlaa 
2U  Saaard  Avaaaa 
Utica,  Naa  York 
Attni  Mr.  Attaaall  T.  IhMitb 
Suparviaoiy  Iwgiiiaar 

1  Bandlz  ATlatlea  Gorporatlas 
Haaaarah  IMratMrlaa  OLriaim 
Aauthfiald,  NthlM 
AttRi  M*.  Ralph  M.  Laraan 

1  Mr.  C.R.  AdM 

Phyaiea  Taahnalagy  Dapartaant 
AaraW^aee  OivtalaB 
Tha  Bealag  Ceapaiy 
Saattla  24,  VaHilngtaB 

1  Bryant  ChttchUg  Orlndar  CaapM 
60  Clinton  Aawoa 
AprlngfiaU.  Tanmit 
\ttni  Nr.  Raald  Om 


CadUlaa  Oica  < 

F.O.  Ban  sa^ 

Uatrelt  $,  Nial  ,  ^ 

Attni  Nr.  -2.  Taylar,  FrajaH  Ingiaaar 


I 


t 


1 


1 


1 


1 


X 


X 


X 


X 


do.  af  Copiti 

X  Or.  Oaway  J.  SandaU 
Mraetor  af  DaaalapWBt 
SaiTlar  Raaaarah  aad  DavalapM*  Oa. 
Carrlar  PaHany 
Syraoaaa,  Aav  lark 

I  Chum  Vaught  CarparatiaB 
P.O.  Baa  590r 
Oallaa,  Taaa 
Altai  Mr.  R.C.  Blaylaek 

Viet  Fraaldwt  (Bnglaaartag) 

I  Chryalor  Corparetlan 
Oafcaaa  oparatlaoa 
P.O.  Baa  7$T 
Datroit  31,  WLeUgan 
Atiai  Nr.  C.v.  Boldar 

1  U.O.  Kaaarehar  4  AaaaaUtaa 
Canattlting  BngUaara 

3FOB  Main  Btraat 
Hoaaton  2,  TaM 

2  fr.  B.W.  Binaln^iaB 
Cr^gante  Enginaarlag  laberatoty 
Natlanal  Burtaa  af  Btandarda 
Bauldar,  Celwado 

1  Curtiaa  >  right  Oonparatlaa 
''right  Aarenauttoal  Dlvialaa 
Oapartiwnt  6332 
Voad  Rldga  Naa  JtraMT 
Attni  w.J.  Aarwr 

Chiaf  Fraiaet  Faginaar 

1  Dayatron  Paalfla 

9320  Llneala  Beularard 
laa  Aiigalaa  US,  Oaliferala 
Attni  dohart  I.  Balth 

spaaUl  Fi«3aat  Ingiaaar 

I  Fard  Xnatniaant  CeapM 
3I*10  Thonaan  AvMaia 
lof«  laland  City  I,  Bair  Xnk 
Attni  Hr.  Jarria 

I  Mr.  AdeU  IfU 
Fard  Kotor  OiiNpaiiy 

Bnglaaariag  aad  Raaaarab  Btaff 
r.O.  Bon  20$} 

Dtarbom,  Nlahliaa 

I  Or.  t.  Nayar,  Jr. 

Ben^totalllea  Saatian 


T 


Oaarbam,  Klehigan 


X 


tlM  OwmM  Oigftirttii 
fISX  8*  WijuXnii  iS 
L«  A«|tXM,  QiXlf. 

mat  iwrf  ftlmr,  Kiwrljtr 

“  -  ^ - -  * 


Otawal  Ai«BlM  Mvitlia 

.  CftiAfi 


F.O.  I 


iMm  U,  Q 

Miilfr.  r.v. 


X 


•«ariai  aaX  TiifcHMBl  0«Mv 
0«Mi^  li^iaMrtiV  UkMMiy 
QwttaX  tlMtrU  8«phv 
1  Rlitar  UU 
SthiaMUAri  8w 

Atiat  0.8.  fm,  mm$af  8 


Mr.  UV.  Wlm 

(taMval  BXMtrle  “lapigr 

Aliamlt  Hawawpy  TuMm  OapartMaA 

950  VMiani  Avwm  U4c.  >-7A 

l^na*  IbsaMfeMMAAa  1 

iMMMh  UbaraiarlM 
OmmX  Moiar  Cwfatiaa 
(kB««X  IMara  TmImOmI  Oaatar 
U  Mila  aai  Mourf  8M 
V«ma«  «***^«i— 

AttAi  Ifr.  8.  ialMA  Hmi 

IlwhiBiaal  Dwalipiiat  Om«  1 


A.C«  Spaik  nuc  Difitloa 
OMMral  Moton  CovfdiatUa 
Wli^wk—  1,  HimmtiM 

AttBi  AXlM  KmaAaaa  X 

1^.  VaXtar  Caraw 
laarfott  OlTlalaa 
OaMTAl  hraalaian  laa. 
mo  IMrUa  ivaotta 

Uttla  PAUa,  Mw  JarMgr  I 

Qraaan  Alremft  fnfinaarini  Cart* 

8au«a«a.  Lai^  XaUaA,  law  f  ark 
Attnt  Hr.  Oavld  H.  Cialgf  Jr. 

Haohifilaal  Oaal^i  Saatlaa 

fejittaarlnj  DavMtaaaA  1 


IwtarwilAaaaX  luaimi  I8iali1i>a  Omf» 
taaatrah  t  ikHwitary 
8aa  Jaa#!  teXtfaaaAa 
AtUii  Dr«  V*l.  UiHlaAa 

I«A«  laplarf  li^.  Mv.  Hirt 
Imr,  Ineeaiataial 
110  taala  Amm  M 
OmuA  lMi4a,  lllahliM 


»p.  L.8.  laiT,  T.r. 
laar'  Icaaa  Divlalaa 
Abba  8aad 
Ilyrta,  Oblo 

Dr.  Calua  0.  Ooaiaal,  0/9-90 
Bldf.  801,  naat  8,  HU  AUa 
Loakhaai  Mlaailaa  A  Spaaa  Oa. 
8.0.  B«  90A 
Sumriala,  Califanda 


Dr.  J.  S.  Auaaan 
Littan  SyataM.  Xna. 

$900  Canega  Avanua 
VeaAXaad  uUa,  Gallfamla 

Mr.  Dan  Noara 
Litton  SyatOM 
9900  Caa^  Aaanua 
HoaaUoA  MlUa.  Callfaanla 


Hr,  A.  I.  Tbaaaa 
Aatra  Dlvlaiaa 
Tha  Nu^aardt  Caaporattaa 
16999  Satleajr  Straat 
?an  luya.  Callfamia 

Nr.  Erandall  NrfclM 
Tlea  PvaalAaBt  CHgt) 

NaDemall  Alramft  Carparatlaa 
Trail  in  8t.  -  St.  Laula 
Nunialpal  Alrpawi 
Dck  916 

St.  Laula  9,  Klaaourl 

Or.  Dana  Siamliaht 
Haehanlaal  Taalnaloiy  Tamtiialwl 
1  Hwbart  Oriva 
UthM,  law  T«k 


HTdranautloa.  Xaearparatad 
PlidaU  Sahool  load 
Haward  Caunty 

taural,  Ntrylaal  1 


•  9  • 


Ifr.  J.  I.  Lawar 

Chlaf,  li^iftoar  InarkUA  Ciagiwiti 
HanagniNdl  Aaro  DlalaiaB 
2600  UdMr  KodA 
Mompella.  Hmtaacda 


Nr.  Carl  F.  Qraaaaar,  Jr. 
Dlraator  of  Naaaarah 
■aw  Haapahlra  Sail  laarlBia, 
Patarbaraagh.  lav  la^ablra 


Jaek  A  Holnta,  A  Dtalalcn  af 
Tha  Slaglar  Carp«**blaB 
1789  tf  Stroat 
SulU  909 

Waahlngtan  6,  D»  C. 


Nro.  AUaa  Ward,  Librarian 

Iordan  Dirlalan  of  Ohitad  Alraraft  Cary* 

Halan  Stroat 

■otwalk.  Cannaaticut  1 

lorthrap  Corparatian 
■aralr  Oirlalan 
1001  bat  Iroabay 
KawtharrM,  Gallfamla 

Attn)  TaalwOeal  Xafanatlan,  9189  1 


Nr.  Fbul  A.  Pitt,  Vlaa  Praaldant 
EnglnoarlAf  A  laaaarah 

Solar  Alroraft  Cifiny 
2800  Faalfle  u^war 
San  Dtago  18,  Gallfamla 

Mr.  H.  0.  line 

Sporrjr  Qyraae^  Ca^iqr 

C-8 

Oroat  laab,  law  Tavk 


Hr.  1.  L.  Slmana 
Nortroniaa 

A  Dlrlalan  of  larthrap  Carp. 
900  bat  Oranfatharpa  Avanua 
Anahaln,  Gallfamla 


Nra.  Plaranea  Turmil 
Infinoorlni  Librarian 
Sparry  Qyi^opa  Canpaagr 
Oraat  laak,  law  Tack 


Nartranlea 

A  DlTlalan  of  larthrop  Carp. 

100  Maraa  Stroat 

Norwood,  Wiooaahuaatta 

Attnt  Nr.  B.  L.  Swalnaan,  Taah.  Aoat. 

ProGlalan  PraAaata  Oapartant  1 

Pratt  A  Whitnajr  Alremft 
Divlalon  of  UAC  -  CAHH. 

P.  0.  Ban  611 
Niddlatovn,  Cannaatlaat 

Attn  I  Librarian  1 


Sundotrand  Arlntien-Danvar 
2480  Waat  70th  Aramia 
Danrar  81,  Colorado 

H.  Jama  V.  Salaaal,  Praaldant 
Turboeraft,  Ine. 
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